NIST-GCR-92-614 


Development of a Video Image-Based 
Methodology for Estimating Large Scale Transient 
Hydrocarbon Smoke Plume Size and Extent 


J.T. Leonard, E. K. Budnick, G. C. Back and S. J. Ganey 


Order 0M0415 


NUST 


United States Department of Commerce 
Technology Administration 
National Institute of Standards and Technology 


NIST-GCR-92-614 


Development of a Video Image-Based 
Methodology for Estimating Large Scale 
Hydrocarbon Smoke Plume Size and Extent 


J. T. Leonard 
Naval Research Laboratory 
Washington, DC 


E. K. Budnick, G. C. Back, S. J: Ganey 
Hughes Associates, Inc. 
Columbia, MD 20375 


August 1992 


‘tT OF 
wes Co 
eg 


Sponsored by: 

U.S. Department of Commerce 

Barbara Hackman Franklin, Secretary 

Technology Administration 

Robert M. White, Under Secretary for Technology 
National Institute of Standards and Technology 
John W. Lyons, Director 


Notice 


This report was prepared for the Building and Fire 

Research Laboratory of the National Institute of 

Standards and Technology under purchase:order number OM0415. 
The statements and conclusions contained in this report 

are those of the authors and do not necessarily reflect 

the views of the National Institute of Standards and 
Technology or the Building and Fire Research Laboratory. 


Li 


Development of a Video Image-based 
Methodology for Estimating Large Scale 
Hydrocarbon Smoke Plume Size and Extent 


J.T. Leonard, E.K. Budnick*, G.G. Back*, and S.J. Ganey* 


Naval Research Laboratory 
Washington, DC 20375 


*Hughes Associates, Inc. 
Columbia, MD 21045 


Prepared for 


National Institute of Standards and Technology 
U.S. Department of Commerce 


June 10, 1992 


Abstract 


Interest in burning crude oil as a means to mitigate large scale oil spills at sea led 
to mid-scale evaluation of relevant crude oil burning characteristics. As part of this effort, 
a computer-based field measurement technique for estimating the size, shape, and extent 
of visible smoke plumes was developed. Of the experiments in which measurements 
were made, Test 7 provided data on plume trajectory for the largest distance from the 
pan. Good agreement was obtained between corresponding 35 mm photographs and 
digitized images used to estimate plume volume and trajectory. In Test 7, at 600 s after 
ignition, the leading edge of the smoke plume had risen to 780 m and traveled 1800 m 
downwind from the pan. The estimated total volume of the smoke plume at 600 s after 
ignition was 3.8 x 10° m*. At that time, the total volume of the smoke plume was 
increasing at a rate of 3.2 x 10° m*/s. Limited evaluation indicates that the MS-DOS based 
method provides reasonably accurate estimates of visible snoke plume geometry in the 
near-field. Field accuracy depends on plume size, wind speed and direction, and the 
resolution of the equipment. 
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Development of a Video Imaging-based 
Methodology for Estimating Large Scale 
Hydrocarbon Smoke Plume Size and Extent 


1.0 INTRODUCTION 
1 General 


Recent oil spills, such as the Exxon Valdez in Alaska, have increased the public’s 
awareness of the damage done to the environment by such accidents. One of the critical 
factors in equating the amount of damage done to the environment has been the spill 
response method used to clean up the oil. In the case of the Exxon Valdez, the main spill 
response method was containment of the oil spill with booms and the use of skimmers 
to recover the oil from the water surface. This method proved futile when weather 
conditions such as high winds and rough sea conditions prevailed. An alternative but 
somewhat controversial practice involves burning the oil spill, thus reducing the threat of 
damage to marine life and contamination of the water. A major limitation to burning oil 
at sea is that as the oil spreads on the water, the sea acts as a very effective heat sink. 
Therefore, it is difficult to raise the temperature of the thin layer of oil high enough to 


achieve ignition. 


The process of burning crude oil on water as a means to mitigate oil spills has 
been investigated and researched by the Building and Fire Research Laboratory (BFRL) 
at the National Institute of Standards and Technology (NIST) for the United States 
Minerals Management Service (MMS), the United States Coast Guard (USCG), and the 


American Petroleum Institute (API) [1-4]. Through the use of small scale experiments and 
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calculations, these studies have increased the understanding of the burning process of 
crude oil on water. These studies have included investigation of the burning rate, fire 
radiation, smoke emission, smoke composition and smoke dispersion in the atmosphere. 
The burning of oil as a spill response method has shown positive results at the small 
scale level by NIST. However, these favorable results must be replicated under large 
scale field conditions in order to gain acceptance as a spill response method by local 


Officials and oil spill response professionals. 


In response to this need, a series of intermediate or mid-scale field tests was 
conducted in Mobile, Alabama at the U.S. Coast Guard Fire and Safety Detachment 
facility. The objective of these tests was to provide data on smoke composition, soot 
yield, smoke concentration, thermal radiation, oil and water temperature, residual oil 
quantities, and smoke plume optical density, volume and position. The results from these 
experiments are being used to quantify the combustion efficiency and the potential 


environmental impact of burning crude oil as a means to mitigate oil spills. 
1.2 Scope and Objectives 

One of the tasks associated with the "mid-scale" field test effort involved 
development of a computer program and related field measurement/imaging techniques 
to analyze individual test records and provide estimates of smoke plume volumes and 
position. This task was assigned to the Naval Research Laboratory (NRL) under contract 
from the National Institute of Standards and Technology (NIST). 


Specific objectives associated with this task included the following: 


(1) | development of a computer based methodology for analyzing images of 


smoke plumes; 


(2) design of a video based field test scheme(s) to obtain plume images; 


(3) evaluation/modification of the field test scheme during the mid-scale field 


tests; and 


(4) preliminary review of the feasibility of adapting commercially available 


atmospheric dispersion models to prediction of snoke plume volume and 


extent. 
1.3. Approach 
1.3.1 General 


Two methods, discrete pattern recognition and spectral signature analysis, 
appeared to have applicability to this problem. Discrete pattern recognition is the most 
widely developed technique, adapting technologies such as gray band discrimination and 
radar signature. The most attractive features associated with this approach include the 
portability and relatively modest cost of field equipment, and the availability of support 
hardware and software for analyzing, storing, and retrieving the test results, and 
manipulation and presentation of the data. However, there are several potential limitations 
associated with this approach. The most important limitation may lie in the accuracy of 
the field measurements. Several sources of error exist due to the lack of uniform plume 
shape, the distances associated with the required field of view, and the difficulty in 
discriminating between clouds and smoke. Another limitation is associated with the 
difficulties of plume edge discrimination, which may necessitate an interactive step on the 
part of the user during analysis of the measurements. This step could be computerized 
at a later time depending on an evaluation of the error introduced in incorporating gray 


band discrimination in the analysis software. 


Spectral signature analysis is under development and has been successfully field 
tested on a limited scale for applications such as sky-cloud discrimination. Typical 
output, however, is limited to the percent cloud cover rather than quantitative measures 
of size, volume, and position. It appears that the use of this technique, based on physical 


optics of light scattering (e.g., Rayleigh and Mie scattering), could provide a more 
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accurate discrimination among sky, clouds, and smoke. However, considerable 
development time would be required to isolate appropriate color separation filters and 
test their ability relative to distinction between optical light effects from clouds and smoke 
plumes of varying densities. In addition, each measurement station would require at least 


three cameras, considerably complicating the field testing. 


1.3.2 Digital Imaging and Analysis 


Based on detailed examination of the two available processing methods, the 
interactive discrete pattern recognition (IDPR) method was selected. While this approach 
has several limitations, the extent of its development and its adaptability to field test 
conditions were considered very attractive. The IDPR method was designed for use with 
color cameras, a PC-based MS-DOS computer system, and an image digitizer. The basic 
system utilizes two video camera inputs in order to capture the cross-section and the 
longitudinal shape of the smoke plume. The methodology is capable of accommodating 
additional camera positions such as overhead in the near-field and downrange beyond 
the near-field mixing and turbulent entrainment zone. However, these locations were not 
examined during the mid-scale test series. Weather conditions and scheduling 


necessitated concentration on two camera positions. 


Several commercially available digitizing packages and related imaging software 
packages were evaluated. The selection of a particular package was based on evaluation 


of MS-DOS compatibility, speed, adaptability, and ease of use. 


1.3.3 Field Measurements 


As discussed above, the basic field measurements consisted of two video 
cameras, positioned as close to 90 degrees from one another as possible. The other 
principle constraint was that the cameras were to be located such that one camera was 
positioned upwind of the smoke plume and the other perpendicular to the smoke plume. 


This preferred orientation of the cameras is illustrated in Fig. 1. 
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Fig. 1 - Preferred camera positions for plume imaging in 
mid-scale crude oil burn tests 


The upwind or longitudinal axis camera position was used to record the width and 
height of the plume. The crosswind or latitudinal axis camera position was used to 
record the width, leading edge, and height of the plume from the other vantage point (see 
Fig. 2). 


Frame accurate time codes were marked on each tape during the tests. Both 
camera views were then synchronized using the corresponding time coded frames. Time 


zero was set at ignition of the crude oil. 


1.3.4 Data Processing 


The video image data were then processed based on pattern recognition through 
differences in intensity, where the plume ideally rises on the back drop of a clear sky. 
The cross-sectional area of the smoke plume at a specified time and position was 
determined by processing the two synchronized video images. These values were then 
used to estimate plume volume, smoke production rate and plume trajectory at each 
specified time step. A computer program was developed and used to analyze and store 


the data in files which can later be retrieved and displayed. 


The field measurement package was not designed to permit data reduction and 
analysis in situ. However, such a capability could be developed if deemed necessary for 


future applications. 


1.3.5 Applicability of Existing Climatic Models 


Computer models have been developed which are used to predict vapor cloud 
and aerosol dispersion down range from atmospheric contamination points. Typical 
applications of such models are focussed at examining dilution effects on gas vapor and 


aerosol concentrations and the ground level deposition rate of contaminants. 
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A preliminary review was performed to identify representative climatic models and 
evaluate the applicability of such models to soot dispersion and deposition from burning 


crude oil smoke plumes. 


2.0 APPLICABILITY OF AVAILABLE COMPUTER CODES 


It is desirable to be able to predict the trajectory, downrange extent, and particulate 
deposition from a large crude oil smoke plume. By necessity, a predictive method must 
account for climatic and geographic effects, the burning characteristics of the fuel (crude 
oil), and the burning rate. The scope of this task supported a cursory review of existing 


computer models in order to determine their applicability. 


A review of the U.S. Environmental Protection Agency (EPA) on-line bulletin board 
and the National Technical Information Service (NTIS) revealed a number of computer- 
based computational models used for prediction of large scale climatic effects. The EPA, 
the U.S. National Oceanographic and Atmospheric Administration (NOAA), and the U.S. 
Department of Energy (DOE) have developed public domain computer models for 
predicting global effects due to pollution, large vapor cloud releases, and climatic 
changes. In addition, several proprietary models are available which are similar in scope 


and capability to the public domain models. 


Accessibility limited the evaluation to the public domain codes. The evaluation 
included review of selected documentation, discussions with current users of the models, 


and an indepth review of one representative model, INPUFF 2.0 [5]. 


In general, the models identified in the review are used for prediction of global or 
large scale effects. The effects are limited primarily to the far-field region, i.e., at locations 
miles downrange of an incident. These pollution or contamination models are 
characterized by a strong dependence on climatic/weather conditions, and output 
generally takes the form of mass concentrations at specified distances downrange from 


a release. 


None of the pollution models reviewed incorporated plume buoyancy effects. 
Discussions with current users confirmed that near-field effects (where buoyancy directly 
affects the height, shape, and size of the smoke plume) are not incorporated in these 
models, and efforts to do so would require significant resources. At points downrange 
where plume buoyancy effects are lost, there is some possibility that models such as 
these could reasonably predict release concentrations. However, the input information 
necessary to run these models necessitates detailed knowledge of the downrange 


location and local conditions. This information is not readily available. 


While a more complete investigation may identify suitable codes, the brief review 
of over thirty models available through the EPA bulletin board and discussions with 
several current users indicate that these models are not suitable for predicting near-field, 
fire buoyancy driven effects. A detailed review of INPUFF, a gaussian based dispersion 


model, supported this conclusion. 


3.0 DESCRIPTION OF EQUIPMENT PACKAGE 


The system consists of two Sony EVO-9100 Hi 8 mm color video cameras for data 
collection; one Sony EVO-9500 Hi 8 mm video playback unit for data viewing and 
retrieval; a Truevision Vidl/O video conversion box to convert S-Video output of playback 
unit to RGB 15 kHz for input to the video digitizer; and one Sony Trinitron color monitor 
used for both tape preview and for digital video interpretation. The digital equipment 
consists of a custom designed 20 mHz Intel 80386 based MS-DOS compatible computer 
with a VGA board and monitor and a Data Translations 2871-HSI Video Digitizing card 


for data capture and processing (see Appendix A for a complete list of equipment). 


4.00 INFORMATION PROCESSING 


4.1. Algorithm 


The first step in analyzing the videotapes recorded during the fire tests involved 
digitizing and manipulating single frames of video. Initially, a scaling factor was 
determined from a known object height (e.g., reference point), which provided a reference 
scale in meters per pixel. This was done for both the upwind camera view and the 
crosswind camera view. Due to CPU time constraints, the plume was digitized every 


fourth pixel. This sped up the editing process and calculation times considerably. 


Once the frame is digitized, the computer program analyzes the image and 
determines which areas represent the smoke plume. This is accomplished by contrasting 
the image intensity between the smoke and the surrounding atmosphere. This step is 


performed automatically by the computer program. 


If the program overestimates the plume size by including clouds and other 
background interference or if it underestimates the plume size, then the user can increase 
or decrease the intensity to edit and modify the plume. Once the user is satisfied with 
the shape of the plume, the program stores the data to a bit map. This is done for both 


views. 


The final step involves the determination of volume, production rate, and trajectory. 
The program performs several geometric calculations, using two different approaches, 
to provide estimates of smoke volume, smoke production rate, and plume trajectory. One 
approach incorporates video images from two camera locations; the alternative method 


incorporates a single-camera image. 


4.2 Input 


The input routine is menu driven and prompts the user for each input element. 


The user is required to input the camera(s) compass headings, distances from the base 


10 


of the fire to the cameras, and the wind direction also in terms of a compass heading 
execution. These data, along with other descriptive information, are stored in a data file 


prior to execution of the program. 


4.3. Calculations 


4.3.1 Two-camera Method 


Ideally, in the two-camera approach, the video cameras are positioned at ninety- 
degree angles to one another as illustrated in Figure 1; one camera upwind from the 
plume and the other camera positioned crosswind (e.g., perpendicular) to the plume. 
The direction of the smoke plume is assumed to be the same as the direction of the wind. 
From each of these two views, the number of horizontal pixels across the smoke plume 
is determined from the digitized frame. These values are multiplied by the scaling factor 
to calculate the widths of the plume from the two camera views, where "a" is equal to the 
width of the plume from the crosswind camera location and "b" is the width determined 
from the upwind camera location (reference Figure 2). The widths of the plume are 


measured at preselected height increments of ten meters. 


In anticipation that configuring the cameras exactly as depicted in Fig. 1 may not 
be possible under field conditions, the program was modified to correct the images for 
other than parallel or right angle orientations. In the event that the cameras cannot be 
positioned parallel (upwind) and perpendicular (crosswind) to the smoke plume, the 
plume width and length are corrected using simple geometric relationships (see Figs. 3 
and 4). 
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Fig. 3 - Geometric correction to estimated plume length 
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Fig. 4 - Geometric correction to estimated plume width 
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The corrected perpendicular or crosswind plume width (a’) is determined using the 


equation: 
= (1) 


corrected smoke plume width (m); 
smoke plume width as viewed through the video camera (m); and 


angle between the crosswind camera view and the smoke plume 


where 


latitudinal axis. 
This correction produces an image more closely representative of a view perpendicular 


to the smoke plume. 


The upwind camera view width (b’) is calculated in the same manner. The upwind 


plume width is determined using the following equation: 
b’ = b (cosa) (2) 


corrected smoke plume width (m); 
smoke plume with as viewed through the video camera (m); and 


angle between the upwind camera view and the smoke plume 


where ekg 


longitudinal axis. 


The corrected widths (a’ and b’) are then substituted into equations (3), (4), and 
(5) to calculate the area, volume, and smoke production rate of the plume. Once the 
extent of the smoke plume is beyond the viewing angle of the camera, the computer 


program uses the last calculated width "a" or "b" in succeeding calculations. 


The cross-sectional areas of the plume are calculated as elliptical sections at each 


height increment using the following equation for the area of an ellipse: 
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Age pean (3) 


4 
where = cross-section area of the smoke plume (m’); 
a= _ plume wicth from the crosswind camera view (perpendicular to the 
plume) (m); and 
b = plume width from the upwind camera view (m). 


The volume of the smoke estimated based on the sum of the cross-sectional areas 
of the plume multiplied by a set height increment of ten meters. The following equation 
is used to calculate the volume: 


(4) 


where n= number of height increments. 


The volume calculation is performed for each frame digitized over a set time step. 
The smoke production rate is then estimated by calculating the difference in the smoke 
volume between two sequential frames divided by the time step. The following equation 


is used to calculate the production rate: 


y = Ue led (5) 
‘ At 
where V = volumetric smoke production rate (m*/min); 


V, = volume of smoke from frame n (m*); 
V 


At = time between frames (min). 


= volume of smoke from frame n-1 (m*); and 


n-1 


As the areas and corresponding volumes are determined for each height 


increment, the sum of the vertical pixels is also recorded by the program. The height of 


15 


the smoke plume is then calculated by multiplying the scale factor by the sum of the 


vertical pixels in the plume. 


4.3.2 Single-camera Method 


The computer program was modified to perform the necessary calculations based 
on images from a single camera. For the one-camera approach, a centerline trajectory 
method is used to perform the calculations. These calculations are performed based on 
video images from the view perpendicular to the direction of the plume (crosswind view). 
After the program has retrieved a previously digitized image, the user is prompted to 
overlay a centerline on the smoke plume. An equation is determined for the line segment 
between each point set by the user. The radius of the plume, which is defined as the 
distance between the centerline and the edge of the plume perpendicular to the 
centerline, is then determined at seven locations (equidistant along the centerline) (see 
Fig. 5). The highest and lowest values are deleted, and the remaining five are averaged 
to produce an average radius for that segment of the plume. Deleting the high and low 
values reduces the problem of inconsistencies in the plume which could skew the overall 
results. The volume of each segment of the plume is determined using a solid cylinder 


approximation. The equation used to calculate the volume of each segment is the 


following: 
V, = s REL, (6) 
where S$ = corresponds to the line segment number; 
V, = volume of cylinder (m°); 


; average radius of the cylinder (m); and 


L, 


The volumes of each segment are then added together over the entire length of the 


length of the cylinder/line segment (m). 


plume to produce the total volume. 
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The centerline method is also used to determine the trajectory of the smoke plume. 
Starting from the base of the fire as the origin, the x and y values along the line segments 
making-up the centerline of the plume are recorded over the entire length of the plume 
(see Fig. 6). These values are then used to represent the trajectory of the plume. The 
maximum y value corresponds to the maximum centerline height of the plume. Again, 


the direction of the smoke plume is assumed to be the same as the wind direction. 


4.4 Output 


Program output includes a summary table, graphical displays for each time step, 
and centerline trajectory coordinates. An example of the summary table is illustrated in 
Table 1 and contains required input data and camera location information. Figure 7 
illustrates typical graphical output which includes the digital image of the plume for both 
camera locations and calculated values for plume height, distance, volume, and smoke 
production rate. Calculated values are provided for both the single and two-camera 
approach. Table 2 illustrates the output form for plume trajectory coordinates. It should 
be remembered that the plume trajectory calculations are based on the images from the 


crosswind camera location only. 
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Test 1 


Test No: 


Date of Test: 
Time of Test: 


Test Description: 


Wind Speed: 


Wind Direction: 


1 
4/16/91 
10:00 


20 ft by 20 ft pool fire 


1.3 m/sec 


170 degrees 


Front View Camera Direction: 90 degrees 
Front View Camera Distance: 120m 
Side View Camera Direction: 0 degrees 
Side View Camera Distance: 120m 
Test Length: 60 seconds 

Time Step: 30 seconds 


Table 1. Typical Field Test Input Data (Test 1) 


Table 2. Tabulated Smoke Plume Distance and Height Trajectory 


Compass Heading 350° 
Elapsed Time = 0:30 (m:s) 


Distance (m) 


0 
0.209612 
0.4192271 
0.6288431 
0.838459 
1.04807 
1.257686 
2.096145 
2.934599 
3.773059 
4.611518 
5.449972 
6.288431 
6.288431 
7.12689 
7.965344 
8.803802 
9.642261 
10.48072 
11.31918 


Output at 30 s Elapsed Time (Test 1) 


Plume Trajectory 
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Height (m) 


0 
6.670783 
13.3498 
20.03703 
26.7325 
33.43619 
40.14809 
47.79309 
55.47512 
63.19416 
70.95023 
78.74332 
86.57344 
86.57344 
91.88812 
97.22748 
102.22748 
107.9803 
113.3937 
118.8318 
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5.0 FIELD TESTS 


5.1. Objectives/General Description 


Mid-scale field tests were conducted at the U.S. Coast Guard facility in Mobile, 
Alabama [6]. These tests were used to determine the optimal approach for viewing and 
analyzing the smoke plume from a range of different camera setups. The first series was 
a set of scoping tests which were conducted by NIST in April 1991. This test series was 
used to test the measurement equipment in the field, to select initial video camera 
positions for viewing the smoke plume, to obtain baseline data for the software algorithm, 
and to prepare and develop a preliminary procedure to be used in the field tests 
scheduled for May 1991. This second series of tests was used to implement the 
procedure developed during the initial testing, to evaluate the algorithm developed for 
computing the smoke plume volume, production rate, and trajectory, and to further 


examine different camera locations. 


A key factor affecting video recording of the smoke plume for these tests was the 
weather. The visibility needed to be clear with no rain. This requirement was necessary 
in order to track the plume and provide a clear video image that could be distinguished 
from background. The positions of the cameras for each test were determined by the 


wind direction. 


The wind direction and speed were the most difficult weather conditions to predict 
and account for in deploying the video cameras. Due to limitations in site availability, it 
was difficult to attain a ninety degree position between the two cameras when the wind 


suddenly shifted direction. 


When the wind speed was low, less than 3 m/s (5 mph), the smoke plume traveled 
nearly straight up. At wind speeds less than 2 m/s, the plume was characterized by 
inverted shapes and pockets of entrained air mixed with the visible smoke. Such effects 


could introduce considerable error into the calculations. At wind speeds greater than 
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around 6 m/s, the plume was rapidly diluted, resulting in difficulty in using the discrete 


pattern recognition method to separate the edges of the plume from the background. 


5.2 General Equipment Locations/Orientation 


The optimum location for the two cameras to document the smoke plume was 
unclear before the test series began. To find the preferred positions, the camera setups 
were varied for each of the tests. These locations are shown in Figure 8. Shifting 
weather patterns and limited site availability complicated location of the camera positions. 
As a result of this and the magnitude of the test fires, it became clear during the test 
series that the cameras should be located as far away from the fire source as possible. 
The loss in accuracy is necessary if the position and extent of the visible smoke plume 
are to be estimated beyond one to two minutes. 


53 Procedure 


The series of mid-scale tests conducted were designed to simulate the burning of 
crude oil on sea water for the purpose of examining instrumentation methods and 
providing data. Video taping of the smoke plume provided data for the development of 


software to calculate plume volume and determine its position. 


Based on the weather forecast for the next day, the Mobile Flight Service Station 
and the Coast Guard ATC were notified that a test is likely. On the day of the test, if 
weather conditions appeared favorable, the test crew moved to Little Sand Island, and 
the airborne weather instrumentation was setup to determine the local weather conditions. 
Subsequently, a decision was made by the test director and safety officer to start the test 
sequence. The time between the start of the test sequence and ignition was set at one 


hour. A detailed description of the test procedure can be found in reference [7]. 


Once the wind speed and direction were known, the video equipment was 


deployed. Preferably, one camera was placed upwind of the plume and the other camera 
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perpendicular to the plume. Once the cameras were setup, the compass headings were 
determined and recorded for the camera positions, and the angle above the horizontal 
was measured. The estimated distance to the fire pan was also recorded for both 
cameras. When the pan was ready for fueling all personnel were informed. Video taping 
commenced at this point so that the ignition of the crude oil would be documented on 
the videotape. The ignition of the oil was used to synchronize the two cameras. The 
cameras remained stationary during the entire test sequence. When the fire was out or 


when the plume was out of the video screen, video recording was terminated. 


5.4 Individual Test Configurations/Setup 


As documented in Table 3, it was difficult to locate the cameras in order to 
replicate the preferred test configuration illustrated in Fig. 1. As previously described, this 
configuration included one camera positioned upwind of the smoke plume providing a 
viewing angle parallel to and centered on the downrange longitudinal axis. It also 
included a second camera positioned perpendicular to the downrange longitudinal axis. 
Ideally, this would result in a 90 degree angle between the two cameras. As already 
discussed, changes in wind direction and limitations in site location limited the ability to 
achieve the preferred configuration. This resulted in the need to provide for adjustments 
or corrections in the digitized images, a capability that was incorporated in the computer 
program (see Section 4.3). As can be seen from the test information in Table 3, 


corrections were required for all eight tests. 


For Test 1, both cameras were located relatively close to the fire. Camera 1 was 
positioned perpendicular to the plume at a compass heading of 90 degrees and camera 
2 was located down wind from the plume at 180 degrees. The wind direction was 350 
degrees at a speed of 1.3 m/s (3.25 mph). After about 30 seconds into the test, the 


smoke plume was at a height beyond the field of view for both camera 1 and camera 2. 


In order to track the plume and maintain it within the video screen’s range for a 


longer period, camera 1 was positioned across the bay on McDuffie Island next to the 
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coal station for Test 2 (see Figure 8). The compass heading for camera 1 was 75 
degrees; camera 2 was positioned at 330 degrees on Little Sand Island. The wind 
direction and speed were 288 degrees and 1-3 m/s (2.5-7.4 mph) respectively. Just prior 
to initiation of Test 2, the wind direction shifted. Due to the lack of site availability, camera 
2 could only be moved along the coal station on such short notice. As a result, camera 
1 and camera 2 were not 90 degrees apart, nor were they perpendicular and parallel to 
the smoke plume. The image of the smoke plume extended beyond the camera range 
for camera 2 after 90 seconds. However, the increased distance to camera 1 provided 


a satisfactory viewing angle for up to four minutes for this test. 


The second series of mid-scale field tests began during the second week in May 
1991. For Test 3, the video cameras were both positioned on McDuffie Island to obtain 
the best view possible, given the varying weather conditions. Camera 1 was placed at 
the center of the coal station with a compass heading of ninety degrees, while camera 
2 was positioned on the south end of the coal station at a compass heading of 30 
degrees (see Figure 8). The wind was 3.3 m/s (8.3 mph) at a compass direction of 340 
degrees. This wind speed was strong enough to bend the plume over. The fire was a 
6.1 x 6.1 m (20 x 20 ft) pan with 900 gallons of crude oil, which allowed a burn time of 
approximately 20 minutes. Camera 2 provided a view of the plume initially traveling 
towards the north. Then, as the plume traveled over the water, it changed to a southerly 
direction down Mobile Bay. After approximately two minutes, the plume image was 
extended beyond the viewing angle for camera 1. However, the image was maintained 


within the viewing angle for camera 2 for the remainder of Test 3. 


In Test 4, camera 1 was located at the center of the coal station and camera 2 was 
located at the south end of the coal station. Camera 1 was at a compass heading of 75 
degrees and camera 2 was at 35 degrees. The wind was at a speed of 2.7 m/s (6.8 
mph) with a compass heading of 320 degrees. Again, due to the shifting wind, it was 
difficult to get the necessary ninety degree angle between the two cameras and the 
appropriate upwind and perpendicular views of the smoke plume. Test 4 used a 6.1 x 
6.1 m (20 x 20 ft) pan fire with 1000 gallons of crude oil. The wind speed was sufficient 


to create a conventional bent-plume shape, but after 4 minutes into the test, the visible 


at 


plume was out of the viewing range for both cameras. After reviewing videotapes of 
these tests, it was decided that the camera locations would be changed again to account 


for the large smoke plume. 


For Test 5, camera 1 was positioned on the northern tip of the coal station island 
and camera 2 was located at the center of the coal station (see Figure 8). Camera 1 had 
a direction of 145 degrees and camera 2 was positioned at a compass heading of 50 
degrees, which gave a 95 degree angle between the two cameras. The wind speed was 
3 m/s (7.5 mph) with a compass heading of 350 degrees. This gave an almost downwind 
view of the plume for camera 1 and an angle of 60 degrees to the smoke plume for 
camera 2. The fire was a 10.7 x 10.7 m (35 x 35 ft) pan fire with 4042 L (1068 gal) of 
crude oil. In the camera 1 view, the plume rose out of the viewing range after two 


minutes; for camera 2, the leading edge of the plume was lost after 2.5 minutes. 


In Test 6, camera 1 was located at the north end of the coal station at a compass 
heading of 145 degrees and camera 2 at the center of the coal station with a direction 
of 42 degrees. A 103 degree angle was formed between the two cameras. The wind 
speed was 10 m/s (25 mph) at a compass direction of 20 degrees. As a result, camera 
2 had an almost upwind view of the smoke plume. Camera 1 was positioned 55 degrees 
to the plume. A 10.7 x 10.7 m (35 x 35 ft) pan fire was used with 3815 L (1008 gal) of 
crude oil burned. Since the wind was gusting very high, the plume immediately began 
to bend and the smoke moved rapidly across the bay. After 2 minutes and 45 seconds, 
the leading edge of the plume was beyond the field of view of camera 2. The time to 
extension of the plume beyond the viewing angle for camera 1 was even shorter at only 


2 minutes. 


Beginning with Test 7, camera 1 was located at the Bay Bridge with a compass 
heading of 250 degrees. This was done in order to maintain the plume inside the viewing 
field of the camera for a longer period of time. Camera 2 was again placed at the center 
of the coal station with a direction of 42 degrees. The two cameras were looking almost 
directly across from one another. The wind speed at test time was 3.0 m/s (7.5 mph) 


with a direction of 350 degrees. For camera 1, there was an almost perpendicular view 
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to the smoke plume while camera 2 had a 52 degree angle with the plume. The size of 
the fire was a 10.7 x 10.7 m (85 x 35 ft) pan with 6562 L (1734 gal) of crude oil, which 
gave a burn time of approximately 20 minutes. The wind was strong enough to bend the 
plume and after only 1 minute and 30 seconds the plume was off the screen for camera 
2. The significant increase in distance to camera 1 resulted in video imaging of the entire 


plume for a longer period during the test. 


Test 8 was the largest fire of the large scale field test series. The pan was 15.2 x 
15.2 m (50 x 50 ft) with 10,598 L (2800 gal) of crude oil, which gave a burn time of 
approximately 18 minutes. The wind speed was 1.5 m/s (3.8 mph), resulting in no 
noticeable bending over of the smoke plume. The wind direction was 290 degrees. 
Camera 1 was located near the Bay Bridge with a compass heading of 240 degrees, and 
camera 2 was positioned at the center of the coal station at 42 degrees. Again, the two 
cameras were essentially facing one another. After 11 minutes into Test 8, the smoke 
plume was out of the viewing range of camera 1; camera 2 lost the plume in less than 


1 minute. 


6.0 FIELD TEST RESULTS 


The video images obtained in each of the eight mid-scale crude oil tests were 
reduced and analyzed. The analysis included frame digitizing of the visible smoke plume 
and subsequent calculation of the plume height, downrange distance, volume, volumetric 
production rate, and centerline position. When appropriate, calculations were based on 


both the single and two-camera methods. 
Typical output from the imaging process is illustrated in Fig. 9. The series of 


images were computed from the output from the crosswind camera position in Test 7 to 


demonstrate the sequential change in height and downrange extent of the visible smoke 
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Fig. 9 - Typical image output from crosswind camera position (Test 7 
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over several time steps. The wind speed associated with this test was measured at 3 m/s 
(7.5 mph). 


Table 4 provides a summary of the plume height, downrange distance, volume, 
and volumetric production rate for both the single and two-camera calculation 
procedures. Coordinates associated with the downrange leading edge of the plume are 
based exclusively on output from the centerline trajectory estimates from the single- 
camera (i.e., crosswind position) method. A graphical representation of the leading edge 
position of the plume at each time step is presented in Fig. 10 for Test 7. The calculated 
values indicate that neutral buoyancy was approached at around 10 minutes after ignition. 
This could not be confirmed based on additional downrange images because the plume 
extended beyond the field of view for the crosswind camera position at a horizontal extent 
of 2000 m (6560 ft). 


Calculations and graphical output are summarized in Appendix B for the eight mid- 
scale tests in which video imaging was attempted. It should be recognized that camera 
locations, wind direction, and fire size varied among tests. The results are an attempt to 
demonstrate the efficacy of such measurements and calculations, and are not intended 


for comparative analysis of the fire or plume characteristics among individual tests. 


7.0 DISCUSSION 
7.1 General 


The methodology developed during the mid-scale crude oil fire tests for smoke 
plume imaging appears promising. The discrete pattern recognition technique provided 
reasonable discrimination between visible smoke and background, generally preserving 
the plume shape and extent. The inclusion of an interactive module in the program 
further enhanced the discrimination capability by permitting the user to adjust the gray 
scales or actually modify the boundary location superimposed on the video image prior 


to image digitizing. 
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The computer program is menu driven and relatively simple to use. It is MS-DOS 
compatible and will produce estimates of plume height, extent, volume, and smoke 
production rate for either the single-camera or two-camera method. Both methods 
employ calculations based on the assumptions that the smoke plumes are elliptical in 
cross-section and relatively uniform along the perimeter boundaries. The error introduced 
can be reduced to some degree through selection of smaller incremental sections. 
However, several other sources of error exist, including the selected scale factor, localized 


temperature and wind, and camera viewing angle relative to the smoke plume geometry. 


While difficult to anticipate prior to the field testing, review of the video images 
indicated that in general, cameras must be located at significant distances from the fire 
source in order to capture the full extent of the plume over several minutes. The 
graphical images provided in Appendix B demonstrate the loss of the smoke plume from 
the field of view of the cameras very quickly when the cameras were positioned within 
1000 m of the fire source. Frequently, this occurred in less than one minute after ignition, 


significantly limiting estimates of smoke plume shape and downrange extent. 


Evaluation of various camera locations and vantage points confirmed that the 
camera distance to the base of the plume corresponds approximately to the downrange 
view of the available camera lens (e.g., 1:1 ratio). This would be expected for the type 
of lens used. Therefore, if one is interested in monitoring the plume for a distance of one 
mile downrange, the crosswind camera must be located approximately one mile away 
from the base of the fire and perpendicular to the prevailing wind direction. The results 
from Test 7 (Appendix B) demonstrate this effect. 


The extreme aspect ratio of the height and downrange extent of the plume and the 
distance from the camera to the base of the fire directly affect the resolution of the video 
image. One means of increasing the resolution would be to adapt relatively sophisticated 
lens technology to minimize the effects of the aspect ratio, therefore reducing the camera 
distance and increasing the relative accuracy of the scale ratio. Another approach would 


be to improve the resolution of the digitizing package. 
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An additional concern associated with the camera distance to downrange viewing 
distance ratio of 1:1 is related to estimates of plume volume and trajectory based on the 
centerline calculation method. The accuracy of the results based on the centerline 
calculation is reduced as the scaling factor is increased (e.g., meters per pixel). And, at 
considerable increases in scaling factor, results may be erroneous. For example, if the 
scale is 100 meters per pixel and the distance is such that the plume is measured at two 
pixels wide, the program cannot obtain a true centerline and will default to either one pixel 
position or the other. This problem can be illustrated by again referring to the images for 
Test 7 (see Appendix B). If the images for the upwind and crosswind camera locations 
are examined during the initial growth of the plume, the relatively small image produced 
by the crosswind camera location results in underestimation of the plume dimensions by 
the centerline method. The two-camera method is also affected in that the resolution of 
the video monitor limits the grid size for the upwind camera location. When the elliptical 
sections are calculated with the two-camera method, the volumes are overestimated due 


to overestimation of the plume width from the upwind camera. 


7.2 Visual Comparison (Overall Size and Shape) 


During Test 7, still photographs were taken at one-minute intervals to serve as a 
comparison to the video film footage recorded perpendicular to the direction of the wind 
(Fig. 11). The videotape was then analyzed using the computer program at the same 
times the still photographs were taken to serve as a side-by-side comparison. The raw 
data (bit maps) were printed for the first five minutes of the test in one-minute intervals 


as shown in Fig. 12. The two figures are shown side-by-side in Fig. 13. 


As shown in Fig. 13, the image the computer/operator has determined to be 
smoke is in reasonable agreement with the still photographs. The only visual difference 
between the two image techniques is a roughness around the edges of the computer 
image which is a function of the resolution of the video digitizing board. The resolution 
of the board views the video image as a 512 x 512 matrix of pixels. During the initial 
stages of the fire, the smoke plume may only fill a small percentage of the screen. As 


previously discussed, when the image is small, the resolution of the system is more likely 
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Elapsed Time = 1:00 Min. 


Elapsed Time = 2:00 Min. 


Elapsed Time = 3:00 Min. 


Elapsed Time = 4:00 Min. 


Elapsed Time = 5:00 Min. 


Fig. 11 - Photographic sequence of smoke plume development 
at one minute intervals (Test 7) 
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Elapsed Time = 1:00 Min. 


Elapsed Time = 2:00 Min. 


Elapsed Time = 3:00 Min. 


Elapsed Time = 4:00 Min. 


Elapsed Time = 5:00 Min. 


Fig. 12 - Digitized plume images at one minute intervals (Test 7) 
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Fig. 13 - Comparison of photographic and digitized video plume images 
at one minute intervals (Test 7) 
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to alter the accuracy of the results. As the smoke volume increases, effects of resolution 
dramatically decrease. 


7.3. Size/Distance Comparison 


Evaluating the relative accuracy of this model was difficult given that no accepted 
reference data were available. However, a comparison was made of the location of the 
smoke plume at a set time into the fire. During Test 7, the wind direction tracked down 
the center of the shipping channel leading into downtown Mobile. Just downwind in the 
shipping lane from Little Sand Island is Mobile Shipyard. Mobile Shipyard contains 
numerous towers and cranes which served as good reference points in estimating the 
actual distance downstream the smoke had traveled. The crane selected as a reference 
location is a huge four-mast assembly used to load "C" containers onto cargo ships. 
Using the navigational charts of Mobile Bay, the crane was determined to be located 
approximately 1250 m from the fire location. The video footage and the still photographs 
both support that the plume was directly above the crane six minutes into Test 7. The 
trajectory data for Test 7 at the six-minute mark (Appendix B) indicates that the leading 
edge of the plume was located 1126 m downwind of the fire. This corresponds to an 
error of roughly 10%. Given the magnitude of contributing factors used to conduct this 


analysis, an error of 10% is viewed as quite acceptable. 


8.0 SUMMARY AND CONCLUSIONS 


The absence of a validated, global scale numerical code that can readily predict 
both near-field and far-field effects for burning crude oil on water enhances the need for 
a set of large scale experimental estimates of selected plume parameters. A set of 
experimental estimates are needed in order to validate newly developed models or 


existing models that have been extensively modified. 


A computer based field scale method was developed and evaluated for obtaining 


such estimates. The two-camera approach appeared to provide the most accurate 
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estimates, but may be difficult to employ due to the sensitivity to camera location. While 
the single-camera approach may sacrifice some accuracy, it is significantly more attractive 


from a field implementation standpoint. 


The image digitizing technique and the computer program perform satisfactorily, 
providing an MS-DOS based, interactive user tool. The program is menu driven and 
permits modifications or corrections to video images produced by discrete pattern 
recognition. It has been made as user-friendly as possible, but some experience in its 
use is still recommended. Program documentation is provided in Appendix C. In 
addition, a detailed user’s manual is provided in Appendix D along with a sample printout 


and summary data file. 


Prior to implementation of this methodology, a series of tests should be performed 
to evaluate accuracy. One approach would be to conduct a series of test with a fixed 
object of known size over a range of distances consistent with downrange distances 
associated with model validation. It should be understood that these distances are limited 
relative to desired accuracy. Such an effort would provide estimates of accuracy related 


to image size and camera distance. 


The following conclusions are drawn based on the limited mid-scale testing and 


evaluation of a smoke plume imaging technique: 


1 visible smoke plumes from large, open burning can be characterized in 


terms of geometry and extent with available imaging technology; 


2. comparisons of digitized images and corresponding 35 mm photographs 
indicate acceptable agreement in qualitative terms of plume shape, size, 


and downrange trajectory; 


3. the lens’ focal characteristics for the cameras selected for this effort (high 
resolution 8 mm) resulted in a 1:1 ratio of camera distance to the fire 


location and downrange field of view. For example, if one desires to 
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monitor the smoke plume geometry for 1 km downrange, the crosswind 


Camera must be positioned at least 1 km from the fire location; 


4. the quality of field measurements is dependent on site availability and 
weather conditions; 


5. the results of this effort are limited to plume images no more than 1.8 km 
(~1 mile) in downrange extent. Longer distances can be evaluated, but at 


additional loss in resolution; and 


6. the methodology is suitable for obtaining a set of data for comparison with 


Output from computational models developed to predict smoke plume 
characteristics. However it is recommended that the accuracy of the 
method be evaluated against large objects of known size beforehand in 


order to determine acceptable field measurement limits. 
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The 


18. 


Equipment List 


. Sony EVO-9100 


Hi 8 mm video camera 


Sony EVO-9100 
Hi 8 mm video camera 


. (8) battery pack NP774 
. (2) AC power adapter AC-V30 
. (2) AV connecting cord (2 phono to 2 phono) 


. (2) connecting cord with connectors 


(4-pin mini-DIN to 4-pin mini-DIN) 


. (2) shoulder mount 


. (2) hand grip 


(2) Ambico - video closeup/wideview lens 


. (2) Silva compass 

. (2) Polycast protractor 

. Sony EVO-9500 Hi 8 mm video playback unit 
. Truevision VID I/O video conversion box 

. Sony Trinitron color monitor 

. 20 mHz Intel 80386-based MS-DOS computer 


. Data Translations 2871-HSI video digitizing card 


Relysis VGA monitor 


Ranging Rangematic 1200 


Serial No. 


202909 


202904 


V0311 


100058 
V1004902 
2655 


2909 
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FLUME TRAJECTURY 
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ELAPSED TIME 00:30 
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3.775859 63.19416 
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6.288431 86.57344 
6.288431 86.57344 
7.12689 91.88812 
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8.803862 162.5915 
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186. 48072 113.3937 
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PLUME TRAJECTORY 


TEST NUMBER 1 
ELAPSED TIME 01:00 


COMPASS HEADING 350 degrees 
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DISTANCE (M) HEIGHT (M) 
G 1) 
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1.257686 40. 14809 
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7.965354 64.81836 
9.642271 71.12989 
11.31918 7749897 
11.37158 84. 12283 
11.42398 94. 74928 
11.47639 103.3783 
11.52879 112.0898 
11.58119 126.644 
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FLUME TRAJECTORY 


TEST NUMBER 2 
ELAPSED TIME 62:00 


COMFASS HEADING 290 dearees 
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39. 38624 148.6678 
66.67255 173.218 
73S.96486 197.7774 
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ELAPSED TIME 06:00 


COMFASS HEADING 290 dearees 


75. 88661 


DISTANCE (‘(M} HEIGHT <M) 
i) a) 
4.167842 8.334257 
8.334884 16.678335 
12.580113 25. 08822 
16.66817 33.34791 
28.83521 41.68943 
25.8822 5@. 832468 
31.25275 62.55184 
37.5833 75.8735 
43.75385 87. 4602801 
50.00441 1@@. 13087 
56.25495 112.6654 
62.5055 125.2842 
62.5855 125.2042 
64.58982 158.7749 
66.67253 176.3484 
68. 75686 201.9247 
78.83956 227.5856 
72.9231 253.0854 


278.6/@1 


Zi VYIWYO Tit YHANYO 


Reescce Peri Wh terre Bansrrrit i errbirrrrer tie 1 Goce coreusies seer eueeeuepueteeeete terete tee eeee cause aus ce er eper ee eset ceees sponse cscs ser mrcssesecess ’ vi r BOO . 
OOTY Weve Wee IVI TN Rannnansnant Satrererny SRARAAAS byrne erry 
Reawenee Ss TEEN TN OT Vv aad a a ans 
Respir nye eyes reeves errr arene aera inate ae eee : 
Ry coreerey ay. Tap REG ADSALDRADSAADSEADRALDIA AAA COUT PIR FUR frank Dasa’ 4 
Pamper nenr nd CaARADDDRARAAARAS TOV OCCT CAAAAAAR BAD AAAAMARABAARAAL 
(aAAAADARRADA rvve aR PMDARBDARARAAAAAS 
PARA AS BA. Cas te Rarrak SARALAAA CAS Cans © g e 
T RAAAAAAAPAARADA SAMA 


Oro anyym ayy yin 


3 
er 
ef 
ae 


un 


3 
R 


avy 
M 
pty a 
pus ta a4 VV’ 
ARRAS SpA aR aS teat a A ay 
PARA DIAL AUDA RAADSIASI IIS RAR LADDER ABI a PARAALARERAARARA RADDA DI (wry) SOR ANTTTTONNIAY COATT ray 
FREE YEEEORYT OTT PERI TETY PFO ES rY P eetTE T 2 


W 68T = SANIT 


GLNNIW Yad SHALGW “ND 6FTSET- PES6ET - GL¥H NOILONGOUd 
SYuaLaW “No T46S08Z BZEGBZLZ > GHOWS JO aWNTON 

Sdalan 46 89F > GONVISIG °X¥W 

Sugaldy b6Z 462 > LHOIGH * X98W 


HOVOUddY ANIT YALNID  HOVOUddY MAIN ONL 


ABA : GWE Gasdl rs 
2: Uda LSdl 


58 


oh 


PLUME Th 44= Cho hey 
TES eNUMEER = 
Pbear oewetiiae Oey 

COMPASS HEADING 232°C gveoreas 

DISTANTE if) ne [BRT 
a) @ 
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TEST NUMBER = 
ELAPSED TIME 1:50 


COMPASS HEALING Bac. dacrees 


DISTANCE (¢f HEIGHT (M: 
4) @ 
16.15215 6.375981 
32. 504351 13.8856 
48.45646 19.88311 
64.69861 27 .81442 
80.76876 34.39753 
96.91292 42. 83244 
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137.2955 56.121235 
177.6737 7ivZ2592Z1 
218.0541 87 .44636 
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FLUME TRAJECTORY 


TESPANEMBER, «oS 
ELAPSED TIME 10:00 


COMPASS HEADING 290 degrees 


DISTANCE (M) HEIGHT (M) 

@ Q 

7 57 

14 113 
32 170 
29 227 
36 284 
43 348 
72 421 
199 501 
129 582 
158 642 
186 743 
215 824 
715 824 
272 932 
330 1040 
387 1149 
444 12598 
532 1348 
559 1477 
768 1519 
968 1544 
1161 i578 
1362 1611 
1563 1646 
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Appendix C 


Program Listings 
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VOLUME p. 1 of 7 


1 DEFINT A,C,0O 
2 DIM AS(80,120) ,AF(60,100) , WF(200) ,WS(200) ,V(2,10) 

3 DIM CROSS(2,2) ,POLY(50,2) ,TR(200,2) 

4 DIM YMIN(50) , YMAX(50) , SLOPE(50) , YINT(50) 

6 DIM CR(500) ,O(500) 

8 SCREEN 9:KEY OFF:CLS:GET (0,0)-(10,10),O0:COLOR 4:LOCATE 1,1:PRINT "*":GET (0,0 
20 QE = 10 

30 P$=0:I%=0:J%=0:DIM ARRAY% (3) 


40 DATA &H55 :' PUSH BP 

50 DATA &HCD, &HO5 :' CDOS INT 5 
60 DATA &HSD :' SDH POP BP 

70 DATA &HCB :' CBH RET FAR 
80 :' 90H NOP 


90 P%=VARPTR(ARRAY%(1)):FOR I%=0 TO 4:READ J%: POKE(P%+I%) ,J¢:NEXT I% 
95 SCREEN 0:CLS: LOCATE 5,5:INPUT" ENTER TEST NUMBER "; NUM 

99 FILES="INDAT"+RIGHTS (STRS$ (1000+NUM) ,3)+".DAT" 

100 ' ***** JNPUT INITIAL INFORMATION ***ee2ex%% 

120 OPEN FILE$ FOR INPUT AS #1 

200 LINE INPUT #1, AS 

201 LINE INPUT #1, BS 

202 LINE INPUT #1, C$ 

203 LINE INPUT #1, DS 

204 LINE INPUT #1, E$ 


300 CLOSE #1 

400 ' ***x* CONVERT TO NUMBERS ****#*% 

410 WSPEED=VAL(D$) 

420 WDIR=VAL(ES) 

430 FCAMDIR=VAL(F$) 

440 FCAMDIS=VAL(GS$) 

450 SCAMDIR=VAL(HS) 

460 SCAMDIS=VAL(IS$) 

470 TESTLEN=VAL(J$) 

480 INC=VAL(KS) 

490 LHARDCOPY=VAL (LS) 

495 SHARDCOPY=VAL (MS) 

500 ' #x*x** DISPLAY INFORMATION *#*xe%x%% 

505 GOSUB 550 

507 SCREEN 9 

510 COLOR 15:LOCATE 1,20:PRINT " VISIBLE SMOKE PLUME CALCULATOR" 

511 COLOR 7:LOCATE 3,5:PRINT " TEST NUMBER : ";NUM 

512 COLOR 7:LOCATE 4,5:PRINT DATE OF TEST : " 

513 COLOR 7:LOCATE 5,5:PRINT " TIME OF TEST: " 

514 COLOR 7:LOCATE 7,5:PRINT " TEST DESCRIPTION : ";C$ 

515 COLOR 7:LOCATE 9,5:PRINT " WIND SPEED : ";DS$;:" m / sec" 

516 COLOR 7:LOCATE 10,5:PRINT " WIND DIRECTION : ";ES$;" degrees" 

517 COLOR 7:LOCATE 12,5:PRINT " FRONT VIEW CAMERA DIRECTION : ";FS$;" degrees" 

518 COLOR 7:LOCATE 13,5:PRINT " FRONT VIEW CAMERA DISTANCE : ";G$;" m" 

" SIDE VIEW CAMERA DIRECTION : ";HS;" degrees" 
" SIDE VIEW CAMERA DISTANCE : ";IS;" m" 

521 COLOR 7:LOCATE 17,5:PRINT " TEST LENGTH ":3$;" seconds" 
1 
Li] 
Lid 


TIME STEP ; meKS; seconds" 
CHECK ALL OF THE ABOVE INFORMATION" 
PRESS C TO CONTINUE OR E TO EXIT PROGRAM" 


526 RS=INKEYS 

527 IF RS$="c" OR RS="C" THEN GOTO 600 
528 IF RS="e" OR RS="E" THEN SYSTEM 

530 GOTO 526 

550 ' *xx*xx* OUTPUT TO A DATA FILE *****% 
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551 
ey: 
553 
554 
555 
556 
557 
558 
559 
560 
561 
562 
563 
564 
565 
566 
567 
568 
600 
700 
800 
1000 
1020 
1030 
1035 
1040 
1045 
1050 
2055 
1060 
1500 
1510 
1511 
1512 
1513 
1514 
1515 
1516 
1517 
1518 
1519 
1520 
1521 
1522 
1524 
1526 
1528 
1000 
1000 
1001 
1002 
1060 
1080 
1100 
1120 
1140 
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FILES="SUM®+RIGHTS (STR$ (1000+NUM) ,3)+".DAT" 
OPEN FILES FOR OUTPUT AS #2 

PRINT #2, 

PRINT #2," VISIBLE SMOKE PLUME CALCULATOR" 

PRINT #2," DATE OF TEST ; 

PRINT #2," TIME OF TEST : ";BS$ 

PRINT #2," TEST DESCRIPTION : ";C$ 

PRINT #2," WIND SPEED : ";D$;" m / sec" 

PRINT #2," WIND DIRECTION : ";E$;" degrees" 

PRINT #2," FRONT VIEW CAMERA DIRECTION : ";FS$;" degrees" 
PRINT #2," FRONT VIEW CAMERA DISTANCE "3GSs""' m" 

PRINT #2," SIDE VIEW CAMERA DIRECTION S;sHS7" degrees" 

PRINT #2," SIDE VIEW CAMERA DISTANCE : ";I$;" m" 

PRINT #2," TEST LENGTH : ";J3$;" seconds" 

PRINT #2," TIME STEP : ";K$;" "seconds" 


IF VAL(E$)>180 THEN N$=STRS (VAL(E$)-180) :GOTO 1000 
NS=STRS (VAL(E$) +180) 

FOR C=INC TO TESTLEN STEP INC 

GOSUB 10000 

GOSUB 50000 

IF LHARDCOPY=1 AND TITLEPAGE=0 THEN GOSUB 1500 

GOSUB 26000 

GOSUB 30000 

NEXT C 

CLOSE #2 

GOTO 95 

6 RHERE PRINTOUT INITDATA kkkkkk&k 

SCREEN 9:CLS:TITLEPAGE=1 

COLOR 15:LOCATE 1,20:PRINT " VISIBLE SMOKE PLUME CALCULATOR" 
COLOR 7:LOCATE 4,5:PRINT ™ DATE OF TEST : ";A$ 

COLOR 7:LOCATE 5,5:PRINT ™ TIME OF TEST : ";BS 

COLOR 7:LOCATE 7,5:PRINT " TEST DESCRIPTION : ";C$ 

COLOR 7:LOCATE 9,5:PRINT ™ WIND SPEED ;: ";DS;" m / sec" 

COLOR 7:LOCATE 10,5:PRINT " WIND DIRECTION : ";ES$;" degrees" 
COLOR 7:LOCATE 12,5:PRINT " FRONT VIEW CAMERA DIRECTION : ";FS;" degrees" 
COLOR 7:LOCATE 13,5:PRINT " FRONT VIEW CAMERA DISTANCE : ";GS;" m" 
COLOR 7:LOCATE 14,5:PRINT ™ SIDE VIEW CAMERA DIRECTION : ";HS;" degrees" 
COLOR 7:LOCATE 15,5:PRINT " SIDE VIEW CAMERA DISTANCE : ";IS$;" nm 
COLOR 7:LOCATE 17,5:PRINT " TEST LENGTH : ":J$;" seconds" 

COLOR 7:LOCATE 18,5:PRINT " TIME STEP : ":KS;" seconds" 

LPRINT CHRS$ (12) 

DEF USRO = VARPTR(ARRAY%(1)): Y=USRO(X) 

RETURN 
O !kkk&e&ee INPUT FRONT v DATA kkekkkakaek 
5 SCREEN 9:CLS:COLOR 4:LOCATE 13,30:PRINT "RETRIEVING DATA" 
O FOR I=1 TO 60:FOR J=1 TO 100:AF(I,J)=0O:NEXT J:NEXT I: TR=0 
O FILEAS = "F" RIGHTS (STRS (1000+NUM) , 3) + RIGHTS (STR$(10000+C), 4) + ".DAT" 
O OPEN FILEAS FOR INPUT AS #1 
O INPUT #1,FSCALE, FCOL, FROW 
O FOR ROW = 1 TO FROW 
0 FOR COL = 1 TO FCOL 
O INPUT #1,AF(COL, ROW) 


11410 W=W+AF (COL, ROW) 
11600 NEXT COL 
11610 WF (ROW) =W*FSCALE: W=0 
11800 NEXT’ ROW 
12000 CLOSE #1 

"RETURN 


13000 
20000 'xekeee INPUT SIDE v DATA #*kkexxxee 


20010 FOR I=1 TO 80:FOR J=1 TO 120:AS(I,J)=0:NEXT J:NEXT I 
20020 FILEA$ = "S" + RIGHTS(STRS$(1000+NUM) ,3) + RIGHTS (STR$(10000+C), 4) + ".DAT 
20130 IF C < 10 THEN FILEAS = "S" + RIGHTS(STRS(C), 1) + ".DAT" 
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20600 
20800 
20900 
21000 
21200 
21400 
21410 
21600 
21610 
21800 
22000 
22010 
23000 
24000 
24010 
24020 
24030 
24040 
24050 
24060 
24070 
25000 
25010 
25020 
25030 
25040 
25050 
25060 
25070 
26000 
26001 
26002 
26004 
26005 
26010 
26100 
26110 
26120 
26200 
26210 
26220 
26300 
26310 
26350 
26400 
26520 
26530 
26540 
26545 
26550 
26560 
26570 
26580 
26590 
27000 
27010 
27100 
27110 
27120 
27200 
27210 
27220 
27300 
27310 
27350 
27400 
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OPEN FILEAS FOR INPUT AS #1 
INPUT #1,SSCALE,SCOL,SROW 
SSCALE=FROW/SROW*FSCALE 
FOR ROW = 1 TO SROW 

FOR COL = 1 TO SCOL 

INPUT #1,AS (COL, ROW) 
W=W+AS (COL, ROW) 

NEXT COL 

WS (ROW) =W*SSCALE : W=0 

NEXT ROW 

CLOSE #1 

CLS:COLOR 15 

RETURN 


‘kazezxk SHOW VIEWS ON SCREEN TO SCALE ****« 

LINE (25,200) -—(275,300) ,3,BF 

FOR ROW = 1 TO FROW 

FOR COL = 1 TO FCOL 

IF AF(COL,ROW)=0 THEN GOTO 24060 

LINE (25+ (COL#2-2) *FSCALE/30-FSCALE/30 , 300-ROW*FSCALE/30-FSCALE/30) -(25+(CO 
NEXT COL 

NEXT ROW 

'ekekk SHOW VIEWS ON SCREEN TO SCALE ***** 

LINE (300,200) -(550,300) ,3,BF 

FOR ROW = 1 TO SROW 

FOR COL = 1 TO SCOL 

IF AS(COL,ROW)=0 THEN GOTO 25060 

Bere eee on Saas oeae ot 2 407 ROWSSSCALEZA OT SOCALE/ oir (A007 


NEXT ROW 
‘keke SHOW VIEWS ON SCREEN TO SCALE **#*# 

IF FSCALE >= SSCALE THEN MD=FSCALE*128 ELSE MD=SSCALE*128 
DFSCALE=FSCALE/ (MD/128) : DSSCALE=SSCALE/ (MD/128) 
IF MD>1000 THEN SPACING=500 ELSE SPACING=100 
CLS:COLOR 15 

LINE (25,172) -(281,300) ,3,BF 

! ##% ADD SCALING LINES **# 

FOR I=SPACING TO MD STEP SPACING 
SF=CINT (I/ (MD/128)) 

' #** X SCALING LINES *** 

'IF SF>FCOL*SCALE THEN GOTO 26300 

LINE (25+SF*2, 300) =(25+SF#2,172) , 15 

! #k& Y SCALING LINES 

'IF SF>FROW*SCALE THEN Goro 26300 

LINE (25, 300-SF) -(281,300-SF) ,15 

NEXT I 

FOR ROW = 1 TO FROW 

FOR COL = 1 TO FCOL 

IF AF(COL,ROW)=0 THEN GOTO 26560 
NCOL=COL*DFSCALE : NROW=ROW* DFSCALE 

LINE (25+NCOL*2-2 , 300-NROW) - (25+NCOL*2, 300-NROW) ,0,BF 
NEXT COL 

NEXT ROW 

LOCATE 23,15:PRINT "CAMERA #1" 

LOCATE 12,12:PRINT "LINES = ";SPACING;" m" 
'#e2ek& SHOW VIEWS ON SCREEN TO SCALE ****s 
LINE (300,172)-(556,300) ,3,BF 

' ##* ADD SCALING LINES *** 

FOR I=SPACING TO MD STEP SPACING 
SF=CINT (I/ (MD/128) ) 

' kk X SCALING LINES *** 

'IF SF>SCOL*SSCALE THEN GOTO 27300 

LINE (300+SF*2, 300) -(300+SF*2,172) ,15 

! #&* Y SCALING LINES *** 

'IF SF>SROW*SSCALE THEN GOTO 27300 

LINE (300, 300-SF) -(556, 300-SF) ,15 

NEXT I 
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27520 FOR ROW = 1 TO SROW 
27530 FOR COL = 1 TO SCOL 

27540 IF AS(COL,ROW)=0 THEN GOTO 27560 

27545 NCOL=COL*DSSCALE: NROW=ROW*DSSCALE 

27550 LINE(300+NCOL*2-2, 300-NROW) -(300+NCOL*2,300-NROW) ,0, BF 

27560 NEXT COL 

27570 NEXT ROW 

27580 LOCATE 23,50:PRINT "CAMERA #2" 

27590 LOCATE 12,46:PRINT "LINES = ";SPACING;" m" 

28500 RETURN 

30000 '! x*aee% VOLUME CALCULATION **#2%22%% 

30001 IF FROW*FSCALE > SROW*SSCALE THEN MH=FROW*FSCALE ELSE MH=SROW*SSCALE 
30005 V=0:DIS=0:MAX=0 

30010 FOR H=0 TO MH STEP 10 

30012 HF=CINT(H/FSCALE) : HS=CINT (H/SSCALE) 

30013 IF WF(HF)=0 AND WS(HS)=0 THEN GOTO 30030 

30014 MAX1=H 

30017 IF WF(HF)=0 AND WS(HS)<>0O THEN WF (HF) =WF (HF-1) 

30018 IF WS(HS)=0 AND WF(HF)<>0 THEN WS(HS)=WS(HS-1) 

30019 IF WF(HF)>DIS1 THEN DIS1=WF (HF) 

30020 IF WS(HS)>DIS1 THEN DIS1=WS (HS) 

30022 V=V+(3.14*WF (HF) *WS (HS) ) *10/4 

30030 NEXT H 

30035 V(1,C\INC)=V 

30037 TS=RIGHTS (STR$(100+C\60) , 2) +8: "+RIGHTS (STR$ (100+ C MOD 60),2) 
30038 LOCATE 1,5:PRINT "TEST NUMBER :NUM 

30040 LOCATE 2,5:PRINT "ELAPSED TIME neTs 

30042 LOCATE 3,20: PRINT "TWO VIEW APPROACH CENTER LINE APPROACH" 
30044 LOCATE 5,5:PRINT "MAX. HEIGHT :" 

30046 LOCATE 5,25:PRINT INT (MAX1) 

30048 LOCATE 5,45:PRINT INT(MAX2) 

30050 LOCATE 5,58:PRINT "METERS" 

30052 LOCATE 6,5:PRINT "MAX. DISTANCE :" 

30054 LOCATE 6,25:PRINT INT(DIS1) 

30056 LOCATE 6,45:PRINT INT(DIS2) 

30058 LOCATE 6,58:PRINT "METERS" 

30060 LOCATE 7,5:PRINT "VOLUME OF SMOKE :" 

30062 LOCATE 7,25:PRINT INT(V(1,C\INC) ) 

30064 LOCATE 7,45:PRINT INT(V(2, C\INC) ) 

30066 LOCATE 7,58:PRINT "CU. METERS" 

30068 LOCATE 8,5:PRINT "PRODUCTION RATE :" 

30070 LOCATE 8,25:PRINT INT((V(1,C\INC) -V(1,C\INC-1) ) *60/INC) 

30072 LOCATE 8,45:PRINT INT((V(2,C\INC)-V(2, C\INC-1) ) *60/INC) 

30074 LOCATE 8,58:PRINT "CU. METERS PER MINUTE" 

30076 PRINT #2, 

30078 PRINT #2, 

30079 PRINT #2,"TEST NUMBER "; 

30080 PRINT #2,"TIME ";TS 

30082 PRINT #2,"TWO VIEW APPROACH CENTER LINE APPROACH" 

30084 PRINT #2,"MAX. HEIGHT : ";INT(MAX1);" *sINT(MAX2);" METERS" 
30086 PRINT #2,"MAX. DISTANCE : ";INT(DIS1) ;" "INT (DIS2) 7" METERS" 
30088 PRINT #2,"VOLUME OF SMOKE : "INT (V(1,C\INC) ); "sINT(V(2,C\INC) );" 
30090 PRINT #2,"PRODUCTION RATE :; " INT ((V(1,C\INC) -V(1, C\INC-1) ) *60/INC) ;" 
30092 PRINT ?2, 

30094 PRINT #2, 

30098 IF LHARDCOPY=1 THEN GOTO 30130 

30100 RS$=INKEYS$ 

30110 IF R$="c" OR RS="C" THEN GOTO 31000 

30120 GOTO 30100 

30130 LPRINT CHR$(12) 

30140 DEF USRO = VARPTR(ARRAY%(1)): Y=USRO(X) 

31000 '***% TRAJECTORY **#*%%% 


31010 CLS 

31020 PRINT * PLUME TRAJECTORY ~ 
31022 PRINT 

31023 PRINT "* TEST NUMBER ";NUM 
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31024 PRINT " ELAPSED TIME ";T$ 

31025 PRINT 

31026 PRINT " COMPASS HEADING ";NS$;" degrees” 

31028 PRINT 

31030 PRINT " DISTANCE (M) HEIGHT (M)" 
31032 PRINT 

31120 PRINT #2, " PLUME TRAJECTORY I 
31122 PRINT #2, 

31123 PRINT #2, 


31124 PRINT #2, " TEST NUMBER ";NUM 
31125 PRINT #2, " ELAPSED TIME ";TS$ 

31126 PRINT #2, " COMPASS HEADING ";N$;" degrees" 

31128 PRINT #2, 

31130 PRINT #2, " DISTANCE (M) HEIGHT (M)" 


31132 PRINT #2, 

31140 FOR I=1 TO TR 

31150 PRINT USING " FFT EF PoE ag 
31152 PRINT #2,USING" $F "2TR« 
31160 NEXT I 

32000 ' x**x*x*x% TRAJECTORY PRINTOUT ***2 

32010 IF LHARDCOPY=0 THEN GOTO 32200 

32030 LPRINT CHRS (12) 


)?;TR(I,2) 
pies TE See) 


He 


32040 LPRINT PLUME TRAJECTORY 4 

32041 LPRINT 

32042 LPRINT " TEST NUMBER ";NUM 

32043 LPRINT " ELAPSED TIME ";T$ 

32044 LPRINT 

32046 LPRINT " COMPASS HEADING ";N$;" degrees" 

32048 LPRINT 

32050 LPRINT " DISTANCE (M) HEIGHT (M)"™ 


32052 LPRINT 
32060 FOR I=1 TO TR 
32070 LPRINT USING " #82 ";TR(I,1):TR(I,2) 
32080 NEXT I 
32170 RS=INKEY$ 
32180 IF RS$="c" OR R$="C" THEN GOTO 32200 
32190 GOTO 32170 
RETURN 


32200 : 
50000 ' xeeeeeeee DRAW kkkkkekeee 


50001 FOR I=1 TO 50:SLOPE(I)=0: YMIN(I) =0: YMAX(I)=0: YINT (I) =0:POLY(I,1)=0:POLZ(I, 
50010 RADMIN=10000:RADMAX=0 : XO=0 : YO=0 

50020 CROSS (1, 1)=40:CROSS (1,2) =240 

50030 CROSS (2,1)=40:CROSS(2,2)=240 

50040 GOSUB 50560 

50042 FOR B = 1 TO 640 STEP 16: LINE (1, 260)-(B, 349), 0, BF: NEXT B 

50044 COLOR 2:LOCATE 22,18:PRINT "CALCULATING THE VOLUME OF THE PLUME" 

50050 FOR N=1 TO POLY 

50060 'PRINT POLY(N,1) ,POLY(N,2) 

50070 NEXT N 

50080 FOR N=1 TO POLY-2 

50090 'PRINT POLY(N,1),POLY(N,2) ,POLY(N+1,1) , POLY (N+1, 2) 

50100 IF (POLY(N+1,1)-POLY(N,1))=0 THEN POLY(N+1,1)=POLY(N,1)+1 

50110 IF (POLY(N+1,2)—-POLY(N,2))=0 THEN POLY(N+1,2)=POLY(N,2)+1 

50120 IF POLY(N,2)<POLY(N+1,2) THEN YMIN(N)=POLY(N,2) ELSE YMIN(N) =POLY(N+1,2) 
50130 IF POLY(N,2)>POLY(N+1,2) THEN YMAX(N)=POLY(N,2) ELSE YMAX(N) =POLY(N+1,2) 
50140 SLOPE(N)=(POLY (N+1, 2) -POLY(N,2) ) / (POLY (N+1, 1) -POLY(N, 1) ) 

50150 YINT(N) =POLY(N,2)-SLOPE(N) *POLY(N, 1) 

50160 'PRINT YMIN(N) , YMAX(N) ,SLOPE(N) , YINT(N) 

50170 NEXT N 

50175 V=0 

50180 FOR N=1 TO POLY-2 

50185 RAD=0:RADSUM=0 

50190 FOR Y=YMAX(N) TO YMIN(N) STEP ((YMIN(N) -YMAX(N) ) /6) 

50195 PSET((Y-YINT(N) ) /SLOPE(N) ,Y) ,2 

50200 X=(Y-YINT(N) ) /SLOPE(N) 

50201 NSLOPE=1/ ((-SLOPE(N) ) *2) 
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50202 YINT=Y-NSLOPE*X 
50203 GOSUB 50530 'TRAGECTORY 

50204 P1X=X: P2X=X: P1Y=Y: P2Y=Y 

50205 GOSUB 50250 ‘RADIUS LINES 

50206 LINE(P1X, P1Y) -(P2X, P2Y) ,2 

50207 GOSUB 50500 

50210 NEXT Y 

50215 GOSUB 50400 

50220 NEXT N 

50225 V(2,C\INC)=V 

50230 RETURN 

50240 ' 

50250 ' ***** RADIUS SUBROUTINE ****#% 

50255 P1=0:P2=0 

50257 IF NSLOPE<-20 OR NSLOPE>20 THEN S=.1 ELSE S=1 
50260 FOR NX=64 TO FCOL*4+68 STEP S 

50270 NY=NX*NSLOPE+YINT! : PSET (NX,NY) ,2 

50275 COL=CINT ( (NX-64) /4) :ROW=CINT ( (259-NY) /2) ':PRINT COL, ROW 
50276 IF COL<1 OR COL>1+FCOL THEN GOTO 50300 

50278 IF ROW<1 OR ROW>1+FROW THEN GOTO 50300 

50280 IF AF(COL,ROW)=1 AND P1=0 THEN GOTO 50320 
50290 IF AF(COL,ROW)=0 AND Pl=1 AND P2=0 THEN GOTO 50350 
50295 IF Pl=1 AND P2=1 THEN GOTO 50310 

50300 NEXT NX 

50310 RETURN 

50320 ' #*&ee% FIRST POINT *#%% 

50330 P1LX=NX:P1Y=NY:P1=1 

50340 GOTO 50300 

50350 ' *&kk&* SECOND POINT **** 

50360 P2X=NX: P2Y=NY: P2=1 

50365 IF P2X<X-10 THEN P2=0 

50370 GOTO 50300 

50380 ' 

50390 ' 

50400 ' %*&kke CENTER LINE VOLUME CALCULATION #*#x« 
50410 ' 

50415 ' 

50420 LDX=ABS( (POLY (N+1,1)-POLY(N, 1) )/SIN (VAL (ES) ~VAL(FS) )) 
50425 DD=((LDX) ~2-(POLY (N+1,1) -POLY(N,1) ) * 

50430 LDY=ABS (POLY (N+1,2)-POLY(N, 2)) # ( (VAL (6) 4DD) /VAL(GS) ) 
50435 L=((LDX~2+(LDY*2) “2)*.5) *FSCALE/4 

50445 RADAVG=((RADSUM-RADMIN-RADMAX) /5) *FSCALE/4 
50450 V=V+3.14*RADAVG*2*L 

50455 'PRINT V 

50460 RETURN 

50500 ' *kk&* RADIUS CALCULATION RUNNING SUM **#x%« 
50501 RDX=ABS ( (P1X-P2X) /SIN(VAL(E$) -VAL(F$) )) 

50502 DD=((RDX) ~2-(P1X-P2X)*2) *.5 

50503 RDY=ABS (P1¥—P2¥) * ( (VAL(G$) +DD) /VAL(GS) ) 

50504 RAD=((RDX*2+(RDY*2) “2) ~.5)/2 

50505 RADSUM=RADSUM+RAD 

50506 IF RAD<RADMIN THEN RADMIN=RAD 

50507 IF RAD>RADMAX THEN RADMAX=RAD 

50508 RETURN 

50530 ' s*k** TRAJECTORY *##*%% 

50531 TR=TRt1 

50532 IF X0=0 THEN XO=X 

50533 IF YO=0 THEN YO=Y 

50534 TRX“ABS ( (X-XO) /SIN (VAL (ES) -VAL(FS) ) ) 

50535 TDD=((TRX) ~2—(X—XO) “2 

50536 TRYCANS (Y-YO) £4 (¥ab( GS) +TDD) /VAL(GS) ) 

50537 TR(TR, 1)=TRX*FSCALE/4 

50538 TR(TR,2) =TRY*FSCALE/2 

50539 DIS2=TR(TR, 1) 

50540 MAX2=TR(TR, 2) 

50541 RETURN - 
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50560 
50570 
50580 
50585 
50590 
50600 
50610 
50620 
50630 
50680 
50720 
50722 
50724 
50726 
50730 
50740 
50750 
50754 
50758 
50760 
50770 
50780 
50790 
50800 
50810 
50820 
50830 
50840 
50850 
50860 
50870 
50880 
50885 
50890 
50900 
51000 
51010 
51020 
51030 
51040 
51050 
51060 
51070 
51080 
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© RRARRRKK DRAW PLUME RAKKKAREK 

SCREEN 9:KEY OFF:CLS 

FINISHED=0 :POLY=0 

GOSUB 51000 

POLY=POLY+1 

GOSUB 50722 

IF FINISHED=1 THEN RETURN 

POLY (POLY, 1) =CROSS(1,1) 

POLY (POLY, 2) =CROSS (1, 2) 

GOTO 50590 

'GOSUB 51000 

LOCATE 20,20:PRINT "DRAW THE CENTER LINE OF THE PLUME" 

LOCATE 21,25:PRINT "PRESS P TO SET A POINT" 

LOCATE 22,15:PRINT "PRESS C WHEN COMPLETE OR R TO REDO CENTER LINE" 

PUT (CROSS(1,1)-15,CROSS(1,2)-15) ,O,AND 

GET (CROSS (2,1) -15,CROSS(2,2)-15) -(CROSS (2,1)+15,CROSS(2,2)+15),0 

PUT (CROSS (2,1) -3,CROSS(2,2)-5),CR 

FOR POL=2 TO POLY-1 

IF POLY>2 THEN LINE(POLY(POL-1,1) ,POLY(POL-1,2))-(POLY(POL,1) , POLY (POL, 2) ) 
NEXT POL 

CROSS (1,1) =CROSS (2,1) 

CROSS (1,2) =CROSS (2,2) 

RS=INKEYS 

IF RS=CHR$(0)+"M" AND CROSS(2,1)<572 THEN CROSS(2,1)=CROSS(1,1)+4:GOTO 507 
IF RS=CHRS$(0)+"K" AND CROSS(2,1)>69 THEN CROSS(2,1)=CROSS(1,1)-4:GOTO 5073 
IF RS=CHRS$(0)+"P" AND CROSS(2,2)<256 THEN CROSS(2,2)=CROSS(1,2)+2:GOTO 507 
IF RS=CHRS(0)+"H" AND CROSS(2,2)>4 THEN CROSS(2,2)=CROSS(1,2)-2:GOTO 50730 
IF RS="6" AND CROSS(2,1)<561 THEN CROSS(2,1)=CROSS(1,1)+16:GOTO 50730 

IF RS="4" AND CROSS(2,1)>80 THEN CROSS(2,1)=CROSS(1,1)-16:GOTO 50730 

IF RS$="2" AND CROSS(2,2)<250 THEN CROSS(2,2)=CROSS(1,2)+8:GOTO 50730 

IF RS="8" AND CROSS(2,2)>9 THEN CROSS(2,2)=CROSS(1,2)-8:GOTO 50730 

IF RS="P" OR RS="p" THEN RETURN 

IF RS="R" OR RS="r" THEN GOTO 50560 

IF RS="C" OR RS="c" THEN FINISHED=1: RETURN 

GOTO 50770 

& theké&k&ane LARGE VIEW kark&axk 

"LINE (64,1) -(578,259) ,7,B 

FOR ROW = 1 TO FROW 

FOR COL = 1 TO FCOL 

IF AF(COL,ROW)=0 THEN GOTO 51060 

LINE (64+COL*4—3 , 259-ROW*2-2) —(64+COL*4, 259-ROW*2) ,8,BF 

NEXT COL 

NEXT ROW 

RETURN 
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DEFINT A-R, T-Z 

DIM CROSS(2, 2) 

DIM W(512), ADJUST(1000), VALUE(512) 
DIM POLY(25, 


DIM D(128, 128) 


1 
2 
3 
4 
5 DIM YMIN(50), YMAX(50), SLOPE(50), YINT(50), LIN(10), EDGE(10) 
6 
10 


DECLARE FUNCTION AUINIT CDECL 
DECLARE FUNCTION AUSYNC CDECL (BYVAL N AS INTEGER) 
DECLARE FUNCTION AUDISP CDECL (BYVAL N AS INTEGER) 
DECLARE FUNCTION AUSETOVL CDECL ALIAS "au_set_ovl_plns"™ (BYVAL N AS INTEGER) 
DECLARE FUNCTION AUEND CDECL ALIAS "au_end" 
DECLARE FUNCTION AUPASS CDECL 
DECLARE FUNCTION AUFREZ CDECL 
DECLARE FUNCTION AUGTPX CDECL (BYVAL BN AS INTEGER, BYVAL RN AS INTEGER, BYVA 
DECLARE FUNCTION AUPTPX CDECL (BYVAL BN AS INTEGER, BYVAL RN AS INTEGER, BYVA 
DECLARE FUNCTION AUSOVP CDECL (SEG N%) 

DECLARE FUNCTION AULINE CDECL (BYVAL PC AS INTEGER,SEG N%,BYVAL MODE AS INTE 
DECLARE FUNCTION AUGPOS CDECL (BYVAL ROW AS INTEGER,BYVAL COL AS INTEGER, BYV 
DECLARE FUNCTION AUBOX CDECL (BYVAL HT AS INTEGER, BYVAL WD AS INTEGER, BYVAL 
‘Buffer 0 = Intensity 

: 1 = Saturation 

4 2 = Hue 

. 3 = Overlay 


DECLARE FUNCTION AUBCLR CDECL (BYVAL N AS INTEGER) 
GOTO 1000 
1 #te&eke&zt CAPTURE FRAME **%4%%% 


03 = 1 "Enable bit plane 2 (Write value = 64) 
04 =1 ‘Enable bit plane 1 (Write valu2 = 


) 
AUSYNC (1) 
AUPASS 
AUSOVP (01) 
AUFREZ 


' kekkk REQUEST BIT STATUS *#e%2% 

' ROW & COL = 0 TO 511 Z IS THE RETURNED VALUE 1-256 
X = AUGTPX(0, ROW, COL, 1, 2Z) 
RETURN 


' akeek%% ADD A POINT ON TV SCREEN *##*%% 
' #kkekkt y = 128 -> crosshair point y = 64 =-> plume point *****x* 
X = AUGTPX(3, ROW, COL, 1, 

Y = Y OR Yi 

X = AUPTPX(3, ROW, COL, 1, Y) 

X = AUPTPX(3, ROW - 1, COL, 1, Y). 


RETURN 
kkkkkk REMOVE A POINT FROM SCREEN ***xx 
kkehtt y = 127 -> crosshair point y = 191 -> plume point 


t 

t 

X = AUGTPX(3, ROW, COL, 1, Y1) 

Y = Y AND Y1 

X = AUPTPX(3, ROW, COL, 1, Y) 

X = AUPTPX(3, ROW - 1, COL, 1, Y) 
RETURN 

{ ekkee&eaed PSET ON TV SCREEN azekke&ke&k 
Y = 128 

X = AUPTPX(3, Y, X, 1, Y) 

RETURN 

' #kkekkz&z LINE ON TV SCREEN ***%2% 
POS1=Y1 

POS2=X1 
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956 
958 
960 


970 
976 
978 
980 


1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1275 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
5000 
5100 
5101 
5400 
5410 
5415 
5420 
5430 
5440 
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X = AUGPOS( Y, X, 3) 
X = AULINE(1, POS1,D) 
RETURN 


' a&kekee BOX ON TV SCREEN **%%2 
X = AUGPOS( Y1, X1, 3) 

X = AUBOX(Y2-Y1,X2-X1,D) 

RETURN 


' gee INPUT INITIAL INFORMATION *eteeeee% 


KEY OFF:SCREEN 9:CLS:LOCATE 5,5:INPUT" ENTER TEST NUMBER "; NUM 


FILES="INDAT"+RIGHTS (STRS$ (2000+NDM) , 3)+". DAT" 
OPEN FILES FOR INPUT AS #1 


LINE INPUT #1, IS 

LINE INPUT #1, JS 

LINE INPUT #1, KS 

CLOSE #1 

! &ke*%k% CONVERT TO NUMBERS *#**%% 
WSPEED=VAL (D$) 

WDIR=VAL (ES) 

FCAMDIR=VAL(F$) 


FCAMDIS=VAL(G$) 

SCAMDIR=VAL (HS) 

SCAMDIS=VAL (I$) 

TESTLEN=VAL(J$) 

INC=VAL(K$) 

' zekkzkx DISPLAY INFORMATION **e2x%% 
SCREEN 9:cls 


COLOR 15:LOCATE 1,20:PRINT “ NIST FIRE PLUME CALCULATOR" 


COLOR 7:LOCATE 3,5:PRINT " TEST NUMBER : ";NUM 
COLOR 7:LOCATE 4,5:PRINT " DATE OF TEST : “;A$ 
COLOR 7:LOCATE 5,5:PRINT " TIME OF TEST : ";BS$ 


COLOR 7:LOCATE 7,5:PRINT " TESTS DESCRIPTION :; 


“CS 


COLOR 7:LOCATE 9,5:PRINT " WIND SPEED : ";DS$;" m / sec" 
COLOR 7:LOCATE 10,5:PRINT " WIND DIRECTION : ";E$;" degrees" 


COLOR 7: LOCATE 12,5:PRINT 
COLOR 7:LOCATE 13,5:PRINT 
COLOR 7: LOCATE 14,5:PRINT 


TEST LENGTH de 


COLOR 7: LOCATE 17,5:PRINT =; 
TIME STEP : KS 


COLOR 7: LOCATE 18,5:PRINT 


" FRONT VIEW CAMERA DIRECTION : ";FS;" degrees" 
" FRONT VIEW CAMERA DISTANCE : ";GS$;" m" 
" SIDE VIEW CAMERA DIRECTION ; 

COLOR 7:LOCATE 15,5:PRINT " SIDE VIEW CAMERA DISTANCE : ";I$;" m" 
M J$:" seconds" 
Li ° 


";HS;" degrees" 


" seconds" 


COLOR 9:LOCATE 22,5:PRINT " CHECK ALL OF Te ABOVE INFORMATION" 


COLOR 9:LOCATE 23,5:PRINT " PRESS C TO CONTINUE OR E TO EXIT PROGRAM" 
COLOR 15 

RS=INKEYS 

IF R$="c" OR RS$="C" THEN GOTO 1460 

IF R$="e" OR RS$="E" THEN SYSTEM 

GOTO 1420 

CLS 

1 RERKE FRONT VIEW ANALYSIS **#%x* 

'TESTLEN = 4: INC = 2 

SCREEN 9: KEY OFF: CLS 


COLOR 7: LOCATE 20, 1: PRINT " INSERT THE VIDEO TAKEN PERPENDICULAR TO THE 
COLOR 9: LOCATE 22, 5: PRINT " PRESS C TO CONTINUE WHEN READY OR R TO RETUR 


COLOR 15 

R$ = INKEYS 

IF R$ = "c" OR R$ = "C" THEN GOTO 5460 
‘IF RS="r" OR RS$="R" THEN GOTO TOP 
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5450 
5460 
5470 
5500 
5502 
5505 
5510 
5515 
5520 
5530 
5540 
5550 
5560 
5580 
5590 
6000 
6101 
6400 
6410 
6415 
6420 
6430 
6440 
6450 
6460 
6470 
6500 
6502 
6505 
6510 
6515 
6520 
6530 
6540 
6550 
6560 
6580 
6590 
7000 
8000 
8001 
8002 
8003 
8004 
8005 
8010 
8012 
8014 
8016 
8018 
8020 
8022 
8025 
8030 
8031 
8032 
8035 
8037 
8040 
8042 
8044 
8046 
8048 
8050 
8052 
8056 
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GOTO 5420 
FOR B = 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), 0, BF: NEXT B 

GOSUB 8000 

FOR FVIEW = INC TO TESTLEN STEP INC 

TS=RIGHTS (STRS$ (100+fview\60) ,2)+":"+RIGHTS (STR$(100+ fview MOD 60) ,2) 

COLOR 7: LOCATE 20, 20: PRINT " FORWARD VIDEO TO TIME = ";t$ 

COLOR 9: LOCATE 22, 5: PRINT " PRESS C TO CONTINUE WHEN READY OR R TO RE-UR 


COLOR 15 

RS = INKEYS 

IF RS = "c" OR RS = "C" THEN GOTO 5560 

IF RS = "ry" OR RS = "R" THEN GOTO 5000 

GOTO 5520 

FILEAS = "F"+ RIGHTS (STR$(1000+NUM) ,3) + RIGHT$(STRS$(10000+FVIEW), 4) + ".D 
GOSUB 8035 

NEXT FVIEW 

‘ kekk& SIDE VIEW ANALYSIS *#*2%% 

SCREEN 9: KEY OFF: CLS 

COLOR 7: LOCATE 20, 3: PRINT ™ INSERT THE VIDEO TAKEN PARALLEL TO THE WIND 

COLOR 9: LOCATE 22, 5: PRINT " PRESS C TO CONTINUE WHEN READY OR R TO R=ZUR 


COLOR 15 
R$ = INKEYS 

IF R$ = "c" OR RS = "C" THEN GOTO 6460 

‘IF R$="r" OR RS="R" THEN GOTO TOP 

GOTO 6420 

FOR B = 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), 0, BF: NEXT B 

GOSUB 8000 

FOR SVIEW = INC TO TESTLEN STEP INC 

TS=RIGHTS (STRS$(100+sview\60) ,2)+":"+RIGHTS (STR$(100+ sview MOD 60) ,2) 

COLOR 7: LOCATE 20, 20: PRINT " FORWARD VIDEO TO TIME = ";t$ 

COLOR 9: LOCATE 22, 5: PRINT " PRESS C TO CONTINUE WHEN READY OR R TO RETUR 
COLOR 15 

R$ = INKEYS 

IF R$ = "c" OR R$ = "C" THEN GOTO 6560 

IF RS = "yr" OR RS = "R" THEN GOTO 6000 

GOTO 6520 

FILEAS = "S" + RIGHTS (STRS$(1000+NUM) ,3) + RIGHTS (STRS$(10000+SVIEW), 4) - ". 
GOSUB 8035 

NEXT SVIEW 


GOTO 1000 
! akekkxkexke*e INITIAL SCREEN GRAP AND SCALING **kex%%0 


FOR B = 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), 0, BF: NEXT B 
COLOR 7: LOCATE 20, 13: PRINT " ADJUST THE HEIGHT OF THE CURSOR TO THE =ASE 
COLOR 9: LOCATE 22, 28: PRINT " PRESS P WHEN IN POSITION" 

GOSUB 9000 ' **#%%* yorigin AND SCALING SUBROUTINE ***x*%%% 

YORIGIN = CROSS(2, 2): XORIGIN = CROSS(2, 1) 

FOR B = 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), 0, BF: NEXT B 

COLOR 7: LOCATE 20, 16: PRINT " ADJUST THE HEIGHT OF THE CURSOR TO A KNCAN 
COLOR 9: LOCATE 22, 28: PRINT " PRESS P WHEN IN POSITION" 

GOSUB 9050 ' *##*#%% yorigin AND SCALING SUBROUTINE *****%% 

KNOWN = CROSS(2, 2) 

FOR B = 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), 0, BF: NEXT B 

COLOR 9: LOCATE 21, 10: INPUT " ENTER THE HEIGHT OF OBJECT IN METERS"; ==IG 
IF (YORIGIN - KNOWN)<1 THEN GOTO 8500 

SCALE = HEIGHT / (YORIGIN - KNOWN) 

FOR B = 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), 0, BF: NEXT B 


RETURN 

« wkkkkkxkk NORMAL SCREEN GRAP *2ee22%% 

FOR B = 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), O, BF: NEXT B 

COLOR 7: LOCATE 20, 11: PRINT ADJUST THE HEIGHT OF THE CURSOR TO THE UPPE 
COLOR 9: LOCATE 22, 28: PRINT “ PRESS P WHEN IN POSITION" 

GOSUB 9000 ' ##%#*%*% yorigin AND SCALING SUBROUTINE ****sxx 

YMAX = INT(CROSS(2, 2)/4)*4: XMIN = INT(CROSS(2, 1)/4) *4 

FOR B = 1 TO 640 STEP 16: LINE | (1, 238)-(B, 349), 0, BF: NEXT B 

COLOR 7: LOCATE 20, 10: PRINT ADJUST THE HEIGHT OF THE CURSOR TO THE LOWE 
COLOR 9: LOCATE 22, 28: PRINT " PRESS P WHEN IN POSITION" 

GOSUB 40000 ' rorre ry. yorigin AND SCALING SUBROUTINE *****s* 


130 


VIDEO p. 4 of 7 


8058 YMIN = INT(CROSS(2, 2)/4)*4: XMAX = INT(CROSS(2, 1)/4) *4 

8060 FOR B = 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), O, BF: NEXT B 

8100 SEN = 90 

8101 LINE (64, 1)-(576, 237), 15, BF 

8102 FOR B = 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), 0, BF: NEXT B 

8105 SEN = SEN + 10: GOSUB 20005' *** CALL VIDEO AND DISPLAY DATA FILE *** 

8110 COLOR 7: LOCATE 19, 11: PRINT ™ ANALYZE THE VIDEO IMAGE OVERLAYED ON THE AD 
8112 COLOR 9: LOCATE 20, 12: PRINT * PRESS I TO INCREASE OR D TO DECREASE THE SE 
8114 COLOR 9: LOCATE 21, 21: PRINT " PRESS C TO USE THE CAPTUED IMAGE" 

8116 COLOR 9: LOCATE 22, 24: PRINT ™ PRESS M TO MODIFY THE IMAGE" 

8120 COLOR 9: LOCATE 23, 22: PRINT " PRESS R TO REDO FROM THE START" 

8125 RS = INKEYS 

8130 IF R$ = "c" OR RS = "C" THEN GOTO 8200 

8135 IF RS = "i" OR RS = "I" THEN GOTO 8102' *** INCREASE INTENSITY *#*# 

8140 IF RS = "dad" OR RS = "D" THEN GOSUB 21000: GOTO 8110' *** DECREASE INTENSITY 
8145 IF RS = "m" OR RS = "M" THEN GOSUB 31000: GOTO 8180! *** MODIFY THE IMAGE * 
8148 IF RS = "r" OR RS = "R" THEN GOTO 8170! *** REDO FROM THE START *#* 

8150 GOTO 8125 

8170 X = AUBCLR(3) 

8174 GOTO 8100 

8180 IF R=1 THEN GOTO 8170 

8200 (kek&k&eed SAVE DATA REKEKERE 

8205 CLS:LINE(1,1)-(639,300) ,3,BF 

8210 COLOR 4:LOCATE 23,31:PRINT "CAPTURING IMAGE" 

8240 OPEN FILEAS FOR OUTPUT AS #1 

8250 PRINT #1,SCALE*4,1+(XMAX-XMIN) \4, 1+ (YMIN-YMAX) \4 

8260 FOR ROW = YMIN TO YMAX STEP -4 

8270 FOR COL = XMIN TO XMAX STEP a 

8280 PRINT #1,D(COL / 4, ROW / 4); 

8290 IF D(COL / 4, ROW y, 4)=1 THEN LINE (64 + COL, ROW / 2 + 2)-(64 + COL + 3, ( 
8300 NEXT COL 

8305 PRINT #1, 

8310 NEXT ROW 

8320 CLOSE #1 

8325 CLS 

8330 RETURN 

8500 ' #kx*xx** POOR SCALING LOGIC *#*eex% 

8510 BEEP 

8520 FOR B = 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), 0, BF: NEXT B 

8530 COLOR 4: LOCATE 20, 21: PRINT " THE SCALING LOGIC IS INCORRECT" 

8540 COLOR 7: LOCATE 21, 5: PRINT " THE BASE OF THE FIRE AND THE KNOWN HEIGHT CA 
8550 COLOR 9: LOCATE 23, 27: PRINT " PRESS C TO CONTINUE" 

8560 R$ = INKEYS 

8570 IF RS = "c" OR RS = "C" THEN GOTO 8000 

8580 GOTO 8560 

9000 ' x*kexk*e yorigin AND SCALING SUBROUTINE ******* 

9010 ‘SCREEN 9:KEY OFF:CLS 

9020 GOSUB 475 ' **** SCREEN GRAB *#**#2#% 

9030 CROSS(1, 1) = 80: CROSS(1, 2) = 246 

9040 CROSS(2, 1) = 80: CROSS(2, 2) = 246 

9050 'LINE(64,1)-(577,258) ,7,B 

9060 'LINE(CROSS(1,1)-10, CROSS (1,2) )—(CROSS(1,1)+10,CROSS (1,2) ) 
9070 'LINE(CROSS(1,1) ,CROSS(1,2)-10)-(CROSS (1,1) , CROSS (1, 2)+10) 
9080 'LINE(CROSS(2,1)-10, CROSS (2,2) )— (CROSS (2,1)+10,CROSS (2,2) ) 
9090 Pena nep 1), CROSS (2, 2)-10)- (CROSS (2,1) , CROSS (2,2)+10) 
9091 Y = 127: = CROSS(1, es ROW = CROSS(1, 2): GOSUB 850 


9092 Y 127: COL = CROSS(1, 1) + 2: ROW = CROSS(1, 2): GOSUB 850 
9093 Y = 127: COL = CROSS(1, 1) - 2: ROW = CROSS(1, 2): GOSUB 850 
9094 Y = 127: COL = CROSS(1, 1): ROW = CROSS(1, 2) + 2: GOSUB 850 
9095 Y = 127: COL = CROSS(1, 1): ROW = CROSS(1, 2) - 2: GOSUB 850 
9096 Y = 128: COL = CROSS(2, 1): ROW = CROSS(2, 2): GOSUB 800 
9097 Y = 128: COL = CROSS(2, 1) + 2: ROW = CROSS(2, 2): GOSUB 800 
9098 Y = 128: COL = CROSS(2, 1) - 2: ROW = CROSS(2, 2): GOSUB 800 
9099 Y = 128: COL = CROSS(2, 1): ROW = CROSS(2, 2) + 2: GOSUB 800 
9100 Y = 128: COL = CROSS(2, 1): ROW = CROSS(2, 2) - 2: GOSUB 800 
9105 CROSS(1, 1) = CROSS(2, 1) 
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9110 CROSS(1, 2) = CROSS(2, 2) 

9120 R$ = INKEYS 

9170 IF R$ = CHRS(0) + "M" AND CROSS(2, 1) < 506 THEN CROSS(2, 1) = CROSS(1, 1) 
9175 IF RS = CHRS(0) + "K* AND CROSS(2, 1) > 5 THEN CROSS(2, 1) = CROSS(1, 1) - 
9180 IF R$ = CHRS(0) + “"P® AND CROSS(2, 2) < 509 THEN CROSS(2, 2) = CROSS(1, 2) 
9185 IF R$ = CHRS$(0) + "H"® AND CROSS(2, 2) > 3 THEN CROSS(2, 2) = CROSS(1, 2) = 
9190 IF R$ = "6" AND CROSS(2, 1) < 495 THEN CROSS(2, 1) = CROSS(1, 1) + 16: GOTO 
9195 IF R$ = "4" AND CROSS(2, 1) > 16 THEN CROSS(2, 1) = CROSS(1, 1) - 16: GOTO 
9200 IF R$ = "2" AND CROSS(2, 2) < 495 THEN CROSS(2, 2) = CROSS(1, 2) + 16: GOTO 
9205 IF R$ = "8" AND CROSS(2, 2) > 16 THEN CROSS(2, 2) = CROSS(1, 2) - 16: GOTO 
9210 IF R$ = "P" OR RS = "p" THEN GOTO 9300 

9220 GOTO 9100 

9300 Y = 127: COL = CROSS(1, 1): ROW = CROSS(1, 2): GOSUB 850 

9302 Y = 127: COL = CROSS(1, 1) + 2: ROW = CROSS(1, 2): GOSUB 850 

9304 Y = 127: COL = CROSS(1, 1) - 2: ROW = CROSS(1, 2): GOSUB 850 

9306 Y = 127: COL = CROSS(1, 1): ROW = CROSS(1, 2) + 2: GOSUB 850 

9308 Y = 127: COL = CROSS(1, 1): ROW = CROSS(1, 2) - 2: GOSUB 850 

9310 RETURN 

20000 ' *** CALL VIDEO AND DISPLAY DATA FILE *** 

20002 GOSUB 475 ' **** SCREEN GRAB *#*x«2% 

20004 LINE (64, 1)- (ie: 257), 15, BF 

20005 FOR I = 1 TO 128 

20010 FOR J = 1 TO 128 

20015 D(i,j)=0 

20020 NEXT J 

20025 NEXT I 

20040 FOR ROW = YMIN TO YMAX STEP -4 

20050 FOR COL = XMIN TO XMAX STEP 4 

20070 GOSUB 700 

20080 IF Z > SEN THEN GOTO 20100 

20085 D(COL / 4, ROW / 4) =1 

20088 Y = 64: GOSUB 800 

20095 LINE (64 + COL, ROW / 2 + 2)-(64 + COL + 3, (ROW + 3) / 2 + 2), 0, BF 
20100 NEXT COL 

20120 NEXT ROW 

20130 RETURN 

21000 '* DECREASE INTENSITY 

21002 SEN = SEN - 10 

21040 FOR ROW = YMIN TO YMAX STEP -4 

21050 FOR COL = XMIN TO XMAX STEP 4 

21070 GOSUB 700 

21080 IF Z > SEN THEN GOTO 21088 ELSE GOTO 21100 

21088 Y = 191: GOSUB 850 

21090 D(COL / 4, ROW / 4) = 0 

21095 LINE (64 + COL, ROW / 2 + 2)-(64 + COL + 3, (ROW + 3) / 2 + 2), 15, BF 
21100 NEXT COL 

21120 NEXT ROW 

21125 

31000 © BRERERE MODIFY IMAGE kekaz&kxe&k&keak 

31001 CROSS(1, 1) = XMIN: CROSS(1, 2) = YMIN 

31002 CROSS(2, 1) = XMIN: CROSS(2, 2) = YMIN 

31005 FOR B = 1 TO 640 STEP 16: LINE (1, 238)-(B, 849) 0, BF: NEXT B 

31010 COLOR 7: LOCATE 20, 18: PRINT " MOVE THE CURSOR TO THE DESIRED LOCATION" 
31015 COLOR 9: LOCATE 21, 12: PRINT " PRESS -B- TO BOX -A= TO ADD OR -E- TO ERAS 
31020 COLOR 9: LOCATE 22, 22: PRINT " PRESS C TO USE THE CAPTUED IMAGE" 

31022 COLOR 9: LOCATE 23, 23: PRINT " PRESS R TO REDO FROM THE START" 

31025 Y = 127: COL = CROSS(1, 1): ROW = CROSS(1, 2): GOSUB 850 

31030 Y = 127: COL = CROSS(1, 1) + 2: ROW = CROSS(1, 2): GOSUB 850 

31035 Y = 127: COL = CROSS(1, 1) - 2: ROW = CROSS(1, 2): GOSUB 3850 

31040 Y = 127: COL = CROSS(1, 1): ROW = CROSS(1, 2) + 2: GOSUB 850 

31045 Y = 127: COL = CROSS(1, 1): ROW = CROSS(1, 2) - 2: GOSUB 850 

31050 Y = 128: COL = CROSS(2, 1): ROW = CROSS(2, 2): GOSUB 800 

31055 Y = 128: COL = CROSS(2, 1) + 2: ROW = CROSS(2, 2): GOSUB 800 

31060 Y = 128: COL = CROSS(2, 1) - 2: ROW = CROSS(2, 2): GOSUB 800 

31065 Y = 128: COL = CROSS(2, 1): ROW = CROSS(2, 2) + 2: GOSUB 800 

31070 Y = 128: COL = CROSS(2, 1): ROW = CROSS(2, 2) - 2: GOSUB 800 
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31075 
31080 
31085 
31090 
31095 
31100 
31105 
31110 
31115 
31120 
31125 
31130 
31132 
31135 
31140 
31145 
31150 
31500 
31505 
31510 
31515 
31520 
31600 
31605 
31610 
31615 
31620 


40000 
40005 
40010 
40020 
40030 
40040 
40220 
40230 
40240 
40250 
40260 
40270 
40280 
40290 
40300 
40310 
40320 
40330 
40340 


41000 
41010 
41020 
41030 
41040 
41050 
41060 
41070 
41080 
41090 
41100 
41110 
41120 
41130 
41140 
41150 
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CROSS(1, 1) = CROSS(2, 1) 
CROSS(1, 2) = CROSS(2, 2) 


R=0: RS=INKEYS ) 

IF R$ = CHRS$(0) + "M" AND CROSS(2, 1) < 506 THEN CROSS(2, 1) = CROSS(1, 1) 
IF RS = CHRS(0) + "K" AND CROSS(2, 1) > 5 THEN CROSS(2, 1) = CROSS(1, 1) - 
IF RS = CHRS(0) + "P" AND CROSS(2, 2) < 509 THEN CROSS(2, 2) = CROSS(1, 2) 
IF RS = CHRS(0) + "H" AND CROSS(2, 2) > 3 THEN CROSS(2, 2) = CROSS(1, 2) - 
IF R$ = "6" AND CROSS(2, 1) < 495 THEN CROSS(2, 1) = CROSS(1, 1) + 16: GOT 
IF R$ = "4" AND CROSS(2, 1) > 16 THEN CROSS(2, 1) = CROSS(1, 1) - 16: GOTO 
IF R$ = "2" AND CROSS(2, 2) < 495 THEN CROSS(2, 2) = CROSS(1, 2) + 16: GOT 
IF R$ = "8" AND CROSS(2, 2) > 16 THEN CROSS(2, 2) = CROSS(1, 2) - 16: GOTO 
IF R$ = "c" OR R$ = "C" THEN RETURN 

IF R$ = "b" OR RS = "B" THEN GOTO 41000! **#* BOX #2 

IF R$ = "a" OR R$ = "A" THEN GOTO 31600' *** ADD A POINT *** 

IF R$ = "e" OR R$ = "E" THEN GOTO 31500' **#* ERASE A POINT *** 

IF R$ = "r" OR R$ = "R" THEN R=1 :RETURN! *** REDO FROM THE START *** 
GOTO 31085 


( wke&eek ERASE A POINT **2z% 

Y = 191: COL = CROSS(1, 1): ROW = CROSS(1, 2): GOSUB 850 
D(COL / 4, ROW / 4) = 0 
LINE (64 + COL, ROW / 2 + 2)-(64 + COL + 3, (ROW + 3) / 2 + 2), 15, BF 
GOTO 31085 

& wkke ADD A POINT kkakk 

Y = 64: COL = CROSS(1, 1): ROW = CROSS(1, 2): GOSUB 800 
D(COL / 4, ROW / 4) = 1 
LINE (64 + COL, ROW / 2 + 2)-(64 + COL + 3, (ROW + 3) / 2 + 2), 0, BF 
GOTO 31085 


© RRHAERERE BOX DRAW SUBROUTINE *2x*xe2% 
X1=XMIN : Y1=YMAX 

D=0: X2=CROSS (1,1) :¥2=CROSS(1,2) :GOSUB 970 
D=1:X2=CROSS (2,1) :¥2=CROSS (2,2) :GOSUB 970 
CROSS(1, 1) = CROSS(2, 1) 

CROSS(1, 2) = CROSS(2, 2) 


RS = INKEYS 

IF R$ = CHRS(0) + "M" AND CROSS(2, 1) < 506 THEN CROSS(2, 1) = CROSS(1, 1) 
IF RS = CHRS(0) + "K" AND CROSS(2, 1) > 5 THEN CROSS(2, 1) = CROSS(1, 1) - 
IF RS = CHRS(0) + "P" AND CROSS(2, 2) < 509 THEN CROSS(2, 2) = CROSS(1, 2) 
IF RS = CHRS(0) + "H" AND CROSS(2, 2) > 3 THEN CROSS(2, 2) = CROSS(1, 2) - 
IF RS = "6" AND CROSS(2, 1) < 495 THEN CROSS(2, 1) = CROSS(1, 1) + 32: GOT 
IF R$ = "4" AND CROSS(2, 1) > 16 THEN CROSS(2, 1) = CROSS(1, 1) - 32: GOTO 
IF RS = "2" AND CROSS(2, 2) < 495 THEN CROSS(2, 2) = CROSS(1, 2) + 32: GOT 
IF RS = "8" AND CROSS(2, 2) > 16 THEN CROSS(2, 2) = CROSS(1, 2) - 32: GOTO 
IF RS = "P" OR R$ = "p" THEN GOTO 40330 

GOTO 40220 


D=0: X2=CROSS (1,1) :Y2=CROSS (1,2) :GOSUB 970 
RETURN 


t tkeekezaeee BOX MODIFY *2eee2n% 


FOR B 


COLOR 


= 1 T0-640 STEP 16: LINE (1, 238)-(B, 349), 0, BF: NEXT B 

7: LOCATE 20, 11: PRINT " ADJUST THE HEIGHT OF THE CURSOR TO THE UPP 
9: LOCATE 22, 28: PRINT " PRESS P WHEN IN POSITION" 

9050 © RRERRERR UPPER LEFT kkeek&kak 

= INT(CROSS(2, 2)/4)%*4: NXMIN = INT(CROSS(2, 1)/4) #4 

= 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), 0, BF: NEXT B 

7: LOCATE 20, 10: PRINT " ADJUST THE HEIGHT OF THE CURSOR TO THE LOW 
9: LOCATE 22, 28: PRINT " PRESS P WHEN IN POSITION" 


X1=NXMIN: Y1=NYMAX:GOSUB 40010 ' *#se%%% TOWER RIGHT ***ee%% 


NYMIN 
FOR B 
COLOR 


= INT(CROSS(2, 2)/4)*4: NXMAX = INT(CROSS(2, 1)/4)*4 
= 1 TO 640 STEP 16: LINE (1, 238)-(B, 349), 0, BF: NEXT B 
9: LOCATE 21, 25: PRINT " PRESS -A- TO ADD OR -E- TO ERASE " 


R=0: R$=INKEY$ 
IF R$ = "a" OR R$ = "A" THEN GOTO 41600' *** ADD A POINT *** 
IF RS = "e" OR R$ = "E" THEN GOTO 41500' *** ERASE A POINT *** 


133 


41160 


41500 
41510 
41520 
41530 
41540 
41550 
41560 
41570 
41580 


41600 
41610 
41620 
41630 
41640 
41650 
41660 
41670 
41680 


GOTO 41130 


VIDEO p. 7 of 7 


' #eeee PRASE A BLOCK #eeeare 
FOR ROW = NYMIN TO NYMAX STEP -4 


FOR COL = 
Y = 191: 
D(COL L 4, 
LINE (64 
NEXT COL 
NEXT ROW 
GOTO 31000 


NXMIN TO NXMAX STEP 4 
GOSUB 850 

ROW / 4) = 

+ COL, ROW / 2 + 2)-(64 + COL + 3, 


| thane ADD A BLOCK xeexe2%%% 


FOR ROW = 
FOR COL 2 


NYMIN TO NYMAX STEP -4 
NXMIN TO NXMAX STEP 4 


Y = 64: GOSUB 800 


D(COL / 4, 
LINE (64 
NEXT COL 
NEXT ROW 
GOTO 31000 


ROW / 4) =1 


+ COL, ROW / 2 + 2)-(64 + COL + 3, 
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(ROW + 3) / 2 + 2), 15, BF 


(ROW + 3) / 2 + 2), 0, BF 
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INTRODUCTION 


Considerable research has been conducted regarding the burning of crude oil at 
sea aS a means to mitigate oil spills. A significant part of this research has been 
conducted by the Building and Fire Research Laboratory (BFRL) at the National Institute 
of Standards and Technology (NIST) for the United States Minerals Management Service 
(MMS), the United States Coast Guard (USCG), and the American Petroleum Institute 
(API). Through the use of small scale experiments and a series of large scale field tests, 
the understanding of the burning process of crude oil on water has been increased. 
These studies included investigation of burning rate, fire radiation, smoke emission, 
smoke composition, and smoke dispersion in the atmosphere. Much of the attention has 
focussed on the development of a computer model which could provide a reasonably 
accurate assessment of the volume of smoke produced and the trajectory of the smoke 
plume resulting from an oil spill fire. This model may also be capable of backing-out the 


quantity of oil burning by the volume of smoke being produced. 


A computer model has been developed to digitize 8 mm video images of smoke 
plumes and provide the user with options to analyze and manipulate the plume data. The 
model requires two recorded video images, one upwind and one crosswind of the plume 
as shown in Figure 1. The crosswind camera records the width, leading edge, and height 
of the smoke plume, while the upwind camera records the width and height of the plume. 
The volume of the smoke plume is calculated using both a one-view and a two-view 


approach. 


For the one-camera approach, a center-line-trajectory method is used to calculate 
the volume from the front camera view. The user is asked to input a center line for the 
smoke plume on the video screen. A radius perpendicular to the center line is 
determined to the edge of the smoke plume. The radii along a discrete line segment are 
averaged together, and then the volume of the plume is determined as solid cylindrical 
shapes. These are then summed up to calculate the total volume. The center line 
method is also the driving logic for determining the trajectory of the smoke plume. In the 


two-camera approach, the cameras are positioned approximately 90 degrees apart, one 
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SMOKE 
PLUME 


CROSS-— WIND 
CAMERA 
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WIND 
DIRECTION 


UPWIND CAMERA POSITION 


Plan View 


Fig. 1 - Preferred video camera orientations 
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upwind from the plume and the other perpendicular. From these two views, the cross- 
sectional areas of the plume are calculated as elliptical sections. If the two cameras are 
not perpendicular and parallel to the smoke plume, the values are adjusted accordingly. 
A set height increment is used to divide the plume into equal sections. The volume is 
then calculated from the widths determined by the crosswind and upwind camera using 
the equation for the area of an ellipse and multiplying by the height increment, i.e. V = 
(xab)/4*dy, where "a" is the width of the plume from the crosswind camera and "b" is the 
width from the upwind camera. The "dy" variable is the predetermined height increment. 
If the smoke plume exceeds the video screen for the upwind view, the previously 
determined width is used in the succeeding calculations of the plume volume and 


production rate through the end of the test run. 
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VIDEO TAPING THE PLUME 


As stated previously, the computer model requires two video tapes, one upwind 
and the other crosswind from the plume. This should position the cameras facing 
approximately 90° apart. If conditions do not permit such positioning, then the cameras 
should be placed as close to these conditions as possible. Deviation from these desired 
positions may decrease the accuracy of the results. Both video cameras should be 
positioned far enough away from the fire to assure that the smoke plume does not leave 
the field of view. The following is a good rule of thumb for positioning the video cameras: 
the width of the view is equivalent to the distance from the fire and the height of view is 
approximately 80% of the distance from the fire. An object of a known height must also 
be visible to each camera. This height will later be used to scale the image during the 
voiume calculations. It is recommended that the height of the object be greater than one- 


half the height of the screen. 


Both cameras should be started at least one minute before ignition. This allows 
for easy synchronization of the video images using the ignition as time zero. If possible, 
the cameras should be mounted on a tripod and should not be moved until the test is 


complete (the fire is out). 


GETTING STARTED—SETUP/WIRING DIAGRAM 
A detailed description of the hardware configurations are shown in the manuals for 


the various pieces of hardware making up the system. A quick reference wire diagram 


has been produced from this information and is shown in Fig. 2. 
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Fig. 2 - Setup/wiring diagram 
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RUNNING THE PROGRAMS 


Initial Information 


The computer model requires a description of the ambient conditions and the camera 
locations during the fire. The file can be produced using any full screen editor. Word 
processing packages are not recommended due to their use of control characters. The 


file should be named as follows: 


FILE NAME "INDAT???.DAT" 
i.e., TEST #1 - INDATQ00.DAT 


lf the test number has fewer than three digits, zeros (@) must be present. 
A description of the variables along with a sample data file is listed below. 
Descriptive variables such as the test description may only contain up to 256 characters. 


Numeric variables are in standard SI units and may contain decimal points but are not 


required. 
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VARIABLE SAMPLE DATA FILE 


TEST NO: 

DATE OF TEST (UP TO 256 CHARACTERS) 

TIME OF TEST (UP TO 256 CHARACTERS) 

TEST DESCRIPTION (UP TO 256 CHARACTERS) 

WIND SPEED - M/SEC 

WIND DIRECTION COMPASS HEADING IN DEGREES 

FRONT VIEW CAMERA DIRECTION COMPASS HEADING IN DEGREES 
FRONT VIEW CAMERA DISTANCE - METERS 

SIDE VIEW CAMERA DIRECTION COMPASS HEADING IN DEGREES 
SIDE VIEW CAMERA DISTANCE - METERS 

TEST LENGTH - SECONDS 

TIME STEP - SECONDS 

FLAG FOR PRINTOUT 1-YES 0-NO 
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2 

4/17/91 
10:00 

2&7 m PAN 
2.2 

108 

75 


650 
330 
120 
60 
30 
1 


Analyzing the Video Image 


The video image manipulation and storage program flow chart is shown in Fig. 3. 
A listing of the program is found in Appendix C of the report, "Development of a Video 
Image-based Methodology for Estimating Large Scale Hydrocarbon Smoke Plume Size 
and Extent." To begin the video imaging program, type "VIDEO," and then press ENTER. 
The following is a list of screens in the order in which they will appear as the program is 


executed. The bold-type print are commands displayed on the computer screen. 


ENTER TEST NUMBER ? 


The test number corresponds to the data file previously made containing all of the 
background information. The test number must be an integer from 1-999. After the 
number has been typed, press ENTER to continue. 
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ENTER TEST NUMBER 
INPUT BACKGROUND 
INFORMATION 


INSERT FRONT VIEW 
VIDEO TAPE 
PERPENDICULAR TO 
THE WIND 


DEFINE BASE 
OF FIRE AND 
KNOWN HEIGHT 


ADVANCE VIDEO TO 
DESIRED TIME 


BOX IN PLUME 


DIGITIZE PLUME 


EDIT PLUME 


STORE PLUME 
TO DATA FILE 


INSERT SIDE VIEW 
VIDEO TAPE 
PERPENDICULAR TO 
THE WIND 


DEFINE BASE 
OF FIRE AND 
KNOWN HEIGHT, 


ADVANCE VIDEO TO 
DESIRED TIME 


BOX IN PLUME 


DIGMZE PLUME 


EDIT PLUME 


STORE PLUME 
TO DATA FILE 


Fig. 3 - Video image manipulation and storage program flow chart 


VISIBLE SMOKE PLUME CALCULATOR 


TEST NUMBER: 8 
DATE OF TEST: 5/31/91 
TIME OF TEST: 9:45 


TEST DESCRIPTION: 50 FT BY 50 FT POOL FIRE 


WIND SPEED: 1 m/sec 
WIND DIRECTION: 110 degrees 


FRONT VIEW CAMERA DIRECTION: 240 degrees 
FRONT VIEW CAMERA DISTANCE: 3500 m 
SIDE VIEW CAMERA DIRECTION: 42 degrees 
SIDE VIEW CAMERA DISTANCE: 650 m 


TEST LENGTH: 60 seconds 
TIME STEP: 30 seconds 


CHECK ALL OF THE ABOVE INFORMATION 
PRESS C TO CONTINUE OR E TO EXIT PROGRAM 


This initial screen displays the background data retrieved from the data file. Check the 
information to assure that it is correct. A capital or a small "C" will allow the user to 


continue. A capital or small "E" will exit the user to the operating system to make 


corrections in the data file. 


INSERT THE VIDEO TAKEN PERPENDICULAR TO THE WIND DIRECTION INTO THE 


PRESS C TO CONTINUE WHEN READY OR R TO RETURN TO THE PREVIOUS MENU 
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This screen serves as a prompt to assure that the proper video tape is in the tape player 
and the player is on and operating properly. Press the VIR button under the input 
selection on the TV monitor to assure the video system is working properly. The tape 
should be at the beginning of the test and include a view of an object of known height 
(i.e. blimp). A capital or small "C" will allow the user to continue with the program. A 


Capital or small "R" will return the user to the previous menu. 


Whenever a cursor appears on the screen, the cursor can be moved using the 
arrow keys on the number pad of the keyboard. When the "NUMLOCK'’ key is turned off, 
the cursor moves in small increments, while when turned on the cursor moves in large 


increments. 


ADJUST THE HEIGHT OF THE CURSOR TO THE BASE OF THE FIRE 
PRESS P WHEN IN POSITION 


This is the first screen in scaling the proceeding video images. Move the cursor to the 
base of the fire using the arrow keys as described above. When in position, press a 
Capital or small "P" to proceed to the next screen. 


ADJUST THE HEIGHT OF THE CURSOR TO A KNOWN HEIGHT 
PRESS P WHEN IN POSITION 


It is recommended that the object be over one-half the height of the screen. This will 
produce more accurate results than an object that is only a small percent of the height 
of the screen. Move the cursor using the arrow keys to a known height (i.e., blimp, 


tower, etc.). When in position, press a capital or small "P" to proceed to the next screen. 
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ENTER THE HEIGHT OF OBJECT IN METERS? 


Enter the height of the object marked by the cursor in the previous screen. The height 
of the object should be input in meters. The height may contain a decimal point but is 


not required. Press ENTER to continue. 


THE SCALING LOGIC IS INCORRECT 
THE BASE OF THE FIRE AND THE KNOWN HEIGHT CAN NOT BE THE SAME POINT 
PRESS C TO CONTINUE 


If this screen appears, there is a problem in the scaling logic. Possible causes for such 
problems are as follows: the base of the fire and the known height falling at the same 
vertical location on the screen; or a zero input on the previous screen for the known 


height. A capital or small "C" will allow the user to repeat the scaling process. 


FORWARD VIDEO TO TIME = 00:30 
PRESS C TO CONTINUE WHEN READY OR R TO RETURN TO THE PREVIOUS MENU 


This screen serves as a prompt to adjust the video to the sample location (desired time). 
The time displayed corresponds to the information present in the background data file. 
The time shown on the video player should correspond to the time of ignition plus the 
time above. The video frame grabber is more efficient when selecting a frame from 
moving video. As the tape plays to the desired time, press a capital or small "C" to 


capture the image. A capital or small "R' will return the user to the previous screen. 
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In most cases, it is desirable to stop the video after the screen has been digitized. 
This may speed up the time required to advance or rewind the video to the next sampling 
location. The acquired image is displayed by pressing the ANALOG/RGB button on the 
TV monitor. If no image is displayed, check the wiring using the wiring diagram found 
in this manual and rerun the program from the beginning. 


ADJUST HEIGHT OF THE CURSOR TO THE UPPER LEFT BOUNDARY 
PRESS P WHEN IN POSITION 


This section of the program allows the user to select a specific section of the screen to 
be analyzed. This reduces the time to process the entire image and reduces the 
presence of other dark objects which may be interpreted as the plume. Move the cursor 


to the desired location using the arrow keys and press "P" when in position. 


ADJUST THE HEIGHT OF THE CURSOR TO THE LOWER RIGHT BOUNDARY 
PRESS P WHEN IN POSITION 


Increase the size of the box to include the entire plume using the arrow keys. The tighter 


the box fits around the plume, the less editing will be required later. Press "P" when in 
position. 


ANALYZE THE VIDEO IMAGE OVERLAYED ON THE ADJACENT MONITOR 
PRESS | TO INCREASE OR D TO DECREASE SENSITIVITY 
PRESS C TO USE THE CAPTURED IMAGE 
PRESS M TO MODIFY THE IMAGE 
PRESS R TO REDO FROM THE START 
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Once the plume has been bounded, the computer will analyze the image in the 
box and make a first cut at determining what portion of the video is plume. The 
computer will draw a series of white dots over the section of the video image on the TV 
screen in which it has determined to be the smoke plume. The sensitivity may be 
increased or decreased by pressing "I" or "D" accordingly. As before, the computer will 
recognize both capital or small letters. It is recommended that the user increase the 
intensity until all of the plume has been selected. Extraneous dots may be removed 
during the modification section of the program. If the user is satisfied with the image, 
pressing the letter "C" will cause the image to be stored to a data file for further 
manipulation. If changes are required, the image can be modified by pressing the letter 
"M4" 


MOVE THE CURSOR TO THE DESIRED LOCATION 
PRESS -B- TO BOX -A- TO ADD OR -E- TO ERASE A POINT 
PRESS C TO USE THE CAPTURED IMAGE 
PRESS R TO REDO FROM THE START 


If only a small area outside the plume has been captured, the data can easily be removed 
by moving the cursor to that location and by pressing the letter "E" to erase. The user 
can add a point to the plume in the same manner by pressing an "A." If a larger area 
exists, the user may want to select (box) an entire area to add or delete. This is 
accomplished by pressing the letter "B." 


ADJUST THE HEIGHT OF THE CURSOR TO THE UPPER LEFT BOUNDARY 
PRESS P WHEN IN POSITION 


The procedure for boxing in an area is the same used to box in the plume. Move the 


cursor to the upper left boundary and press "P" when in position. 
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ADJUST THE HEIGHT OF THE CURSOR TO THE LOWER RIGHT BOUNDARY 
PRESS P WHEN IN POSITION 


Extend the box over the desired area and again press "P" when in position. 


PRESS -A- TO ADD OR -E- TO ERASE 


Once the area is bounded, the letter "A" will add the area to the plume or an "E" will 
remove it from the plume. This process may be repeated until all modifications are 
complete. Once the user is satisfied with the image, pressing the letter "C" will cause the 


image to be stored to a data file for further manipulation. 


Once the image has been stored to a file, the computer will prompt the user to advance 
the video to the next time step, and the time capture process is repeated. When all the 
video images required from the first tape are complete. The computer will prompt the 


user to change tapes, and the process is again repeated. 


INSERT THE VIDEO TAKEN PARALLEL TO THE WIND DIRECTION INTO THE PLAYER 
PRESS C TO CONTINUE WHEN READY OR R TO RETURN TO THE PREVIOUS MENU 


A sample Bit Map (data file) developed during this process is shown in Fig. 4. Due to the 
manner in which the matrix of data is handled, the image is stored upside down and 
backwards. The nomenclature for the file is similar to the initial data file. The first 
character stands for front or side, the next three characters are the test number, and the 
final four are the time (seconds) into the fire from which the image was recorded. The 


first number in the data file is the scaling factor. Each number making up the plume 
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# Cols # Rows 


M/pix 


g 


8 


< 
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SBSseGeseGleGesesseseSGesSesSeSGeGeSseSeSeSseSGeGeGseSesGeGaeSsHGeGesescsesGecseseeseecseseseaess 
SB8eGece8ascecdoeeseGewnsGBOoGSeGoGeGevesveaGeeGrazannHGoGesseGeGesieGaeseseeGeaeeeGaeosoeoeseessa 
BSBSSSSSSSGSSSSSGHrtantntnntnnrantnntGSSSSGnttanntantnrtanrnsBGBegegesveGseGeSGsesoaosesesegaesssssessa 
SSGOS SSS SHG HH HAH HHH HAM H HMMM HN HM KHANH MH HMM NHN K KH HANH SSSeeseGSeseeeaeeGacaeecasvaecass 
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E-Ten ok ek wok wok Bok Boek ok Boe eo Eo eee eee ee ee er ee ee ee ee Oe er ee ee eee eee Se 
BSOS OSS SHH HAHAHA NH HN HH MN MAMAN BA HH HH NANNHNHNH MANA MN HM HM NM NHN NH HMMM MMMM KH HMMA OSSeses 
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SBSSSGSOSGSSSESSSHSOSSHSSHF AOS SHA HAHA HHA HH HAHAHAHAHA BH HHA HH HH NBNHHHAHH MH HH KRH MMM HHH RAM H NHN ASSses 
SGSGOSHSHSSSSESHSSSSSHSESSHSSOKDGSM HAHAHAHAHAHA DAHA HHHHHHHANHN NHN HHMNNHMHAHMHMNMM HHMI KH AOSSsa 
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(Seconds) 


F0080030.DAT 
Test # Time 


Fig. 4 - Sample bit map (data file) 
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represents a square with a side of many meters. The next two numbers are the number 
of columns and rows making up the image. The remaining numbers make up the image. 


The area representing the plume is depicted by the "1's in the file. 


Calculating the Volume and Trajectory of the Plume 


The volume and trajectory calculation and output program flow chart is shown in 
Fig. 5. A listing of the program is found in Appendix C of the report, "Development of a 
Video Image-based Methodology for Estimating Large Scale Hydrocarbon Smoke Plume 
Size and Extent." To begin the computation program, type "VOLUME" then press ENTER. 


ENTER TEST NUMBER ? 


As before, the test number corresponds to the data file containing all of the background 
information. Now the stored images of the plume produced using the video imaging 
program should also be present. The test number must be an integer from 1-999. After 


the test number has been typed, press ENTER to continue. 
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START 


ENTER TEST NUMBER 
INPUT BACKGROUND 
INFORMATION 


RETRIEVE BIT MAPS 
DEVELOPED IN 


THE VIDEO 
PROGRAM 


DRAW THE 
@ OF THE PLUME 
ALONG THE 
FRONT 
VIEW 


CALCULATE THE RADIUS 
ALONG THE G, 
THE TRAJECTORY, 
AND THE VOLUME 

OF THE PLUME 


CALCULATE THE VOLUME 
OF THE PLUME 
USING A TWO 

VIEW APPROACH 


DISPLAY 
RESULTS 
(VOLUMES) 


DISPLAY 
TRAJECTORY 


STORE DATA TO 


DATA FILE 


Fig. 5 - Calculation and output program flow chart 
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VISIBLE SMOKE PLUME CALCULATOR 


TEST NO: 8 
DATE OF TEST: 5/31/91 
TIME OF TEST: 9:45 


TEST DESCRIPTION: 50 FT BY 50 FT POOL FIRE 


WIND SPEED: 1 m/sec 
WIND DIRECTION: 110 degrees 


FRONT VIEW CAMERA DIRECTION: 240 degrees 
FRONT VIEW CAMERA DISTANCE: 3500 m 
SIDE VIEW CAMERA DIRECTION: 42 degrees 
SIDE VIEW CAMERA DISTANCE: 650 m 


TEST LENGTH: 60 seconds 
TIME STEP: 30 seconds 


CHECK ALL OF THE ABOVE INFORMATION 
PRESS C TO CONTINUE OR E TO EXIT PROGRAM 


The initial screen displays the background data retrieved from the data file. Check the 
information to assure that it is correct. A capital or a small "C" will allow the user to 
continue. A capital or small "E" will exit the user to the operating system to make 
corrections in the data file. 


RETRIEVING DATA 
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Once the test conditions have been checked, the computer will retrieve the bit maps of 
the two views at the first time step developed using the video program. Once all of the 
data have been retrieved, the user must first draw the center line of the plume on the 


video recorded perpendicular to the direction of the wind. 


DRAW THE CENTER LINE OF THE PLUME 
PRESS P TO SET A POINT 
PRESS C WHEN COMPLETE OR R TO REDO CENTER LINE 


Move the cursor using the arrow keys to draw the center line of the plume starting at the 
base of the fire. Press the letter "P" to set points along the center line. The computer will 
draw the center line as the user selects the points as shown in Fig. 6. If the center line 
does not look correct, press "R' to redraw the center line from the beginning. Press "C" 


when the center line is complete. 


457, 


Fig. 6 - Centerline of the plume 
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Once the center line is complete, the computer will calculate the volume of the plume 
using a one-view approach. If the front view camera is not exactly perpendicular to the 
direction of the wind (plume), the image is first adjusted accordingly. The computer will 
determine the radius of the smoke plume at seven locations along each line segment 
making up the center line as shown in Fig. 7. The high and low values are neglected, 
and the volume of the plume at that location along the center line will be estimated by a 
cylinder of a length equal to the line segment length and a radius of the average of the 
remaining five radii. These cylindrical volumes are then added to produce the volume of 
the plume. The computer will display each radius/diameter as the values are calculated. 


The center line of the plume will also be used to calculate the trajectory of the plume. 


After the single view volume and trajectory calculations are complete, the volume 
of the plume is then calculated using a two-view approach. If the two cameras are not 
at a ninety-degree angle between one another, the two values are adjusted perpendicular 
and parallel to the smoke plume. A set height increment (10 meters) is used to divide 
the plume up into equal sections. The volume is then calculated from the widths 
determined by the front and side cameras using the equation for the area of an ellipse 
and multiplying by the height increment, i.e. V = (zab)/4*dy, where "a" is the width of the 
plume from the front camera and "b" is the width from the side camera (see Fig. 7). The 
"dy" variable is the predetermined height increment. If the smoke plume exceeds the 
video screen for side view, the previously determined width is used in succeeding 


calculations of the plume volume and production rate through the end of the test run. 


Once these two calculation techniques are complete, the results are displayed as 
shown in Fig. 8. The maximum height and distance may vary dramatically between two 
approaches. The reason for this is that the two-view approach uses the edges of the 
plume to determine this information while the one-view uses the center line of the plume. 


The user can advance to the next screen of data by entering the letter "C." 


159 


Fig. 7 - Logic used to calculate the volume in the one-view approach 
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TEST NUMBER : 1 
ELAPSED TIME : 66:38 


TWO VIEW APPROACH 


MAX. HEIGHT : 
MAX. DISTANCE : 
VOLUME OF SMOKE : 
PRODUCTION RATE : 


LINES = 166 m 


ey 


PE re reece ee TTR ORTON TO LY TOOT TY OOTY TENT E TLE TUTTO NUTS UVTI? 


wey 


CAMERA #1 


CENTER LINE APPROACH 


METERS 
_ METERS 
CU. METERS 
CU. METERS PER 


LINES = 186 mn 


tae oe 


~ 


CA RS 
Wh et 


ET I ON 


A TTD 


RIN I 8 WT 


STN 


CAMERA #2 


Fig. 8 - Sample output screen 
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MINUTE 


PLUME TRAJECTORY 


COMPASS HEADING 110 degrees 


DISTANCE M HEIGHT M 
0 0 
0 11 
0 21 
1 32 
1 43 
1 93 
1 64 
2 80 
5 96 
8 112 

11 128 
14 144 
17 160 
17 160 
21 175 
24 190 
28 205 
32 219 
35 234 


This screen is a list of distances and corresponding heights which make up the 
trajectory of the plume. For these calculations, the model assumes the smoke is moving 
in the direction of the wind. The user can proceed with the program by typing in the 
letter "C." At this point, the process is then repeated using the bit maps produced from 
the next step into the fire. A summary data file is produced each time the program is run. 
The file is stored in ASCIl format under the name SUM???.DAT (i.e. SUMQ@1.DAT for Test 


1). A sample printout and summary data file is found in the Appendix. 
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SUGGESTIONS FOR INCREASED ACCURACY 


° Try to capture the entire plume in both video images. 
: Do not move the camera during the test. 
2 lf conditions permit, scale the image with a height of an object over one-half the 


height of the screen. 


‘ If the wind speed is very low (less than 1 m/sec), use angles parallel and 
perpendicular to the wind in the background data file for the location of the 
cameras. 


. If one view is inferior due to problems with scaling or loss of the plume off the 
screen, rearrange the parameters in the initial data file to use the superior image 
as the front view/perpendicular image. This allows the center line volume 


calculation to be conducted using the superior view. 


‘ If the wind speed is very low (less than 1 m/sec) and one view is inferior, use the 


good view for both video images. 
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Appendix 


Sample Printout and Summary Data File 
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FLUME TRAJECTORY 


TEST NUMBER 3S 
elntiseorni ME O03 30 


COMPASS HEADING 2506 degrees 


DISTANCE (i) HEIGHT (M™) 
a 1) 
9.936396 8.889222 
19.87281 17.9468 
29. 8092 27217275 
39. 74561 36.56785 
49. 68261 46.12972 
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PLUME TRAJECTORY 


TES; NUPSER® S 


SWepoe et LM, ,O13.00 
COMPASS HEADING 350 dearees 


DISTANCE (MM: HEIGHT (Mm) 


@ 1) 
11.17844 8.899744 
22. 35688 17. 98889 
33.535355 27.2.26745 
44.71379 36. 73541 
Joe 89225 46.39278 
67.878635 536. 23949 
91.91166 67.9859 
116.7527 78.35321 
141.5937 90. 84142 
166. 4347 102.1585 
191.2757 114. 46866 
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PLUHE TRAJECTORY 


Pe SeNCnioe ir a oo 
ELON O1:50 


COMPASS HEADING 5G dearses 
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DISTANCE (™) HEIGHT (M> 
a @ 
9.93641 72619325 
19.87276 15.38296 
29.88917 25.2989 
29. 74558 31.34315 
49.681935 39.53971 
09.61834 47.88858 
99.6184 47.88058 
109.2004 69.58726 
158.9824 92.817356 
208.6644 117.8115 
298.3464 144.489 
308.8285 172.8581 
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FLUME TRAJECTORY 


TEST NUMBER & 
BEA sey 7 IMEC O2s 00 


COMFASS HEADING 350 dearees 
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DISTANCE (MM) HEIGHT (M) 
@ a 
9.93641 1@.15912 
19. 87282 28.51865 
29. 886923 31.8546 
39. 74564 41.79095 
49. 68285 92.71972 
99.61834 63.840678 
119.2567 81.3259 
178.8551 100.2541 
238.4735 120.6254 
298.0919 142.4398 
357.7183 165.6973 
417.3284 198.3978 
442.1694 211.5748 
467.6104 235.3532 
491.8514 295. 7328 
316.6924 278.7136 
941.5333 362.2958 


VISIBLE SMOKE PLUME CALCULATOR 

TEST NUMBER : 5 

DATE OF TEST : 5/28/91 

TIME OF TEST : 10:00 

TESTS DESCRIPTION : 35 FT BY 35 FT POOL FIRE 
WIND SPEED : 3.0m / sec 

WIND DIRECTION : 0 degrees 

FRONT VIEW CAMERA DIRECTION : 100 degrees 
FRONT VIEW CAMERA DISTANCE : 950 mn 

SIDE VIEW CAMERA DIRECTION : 40 degrees 
SIDE VIEW CAMERA DISTANCE : 650 mn 

TEST LENGTH ; 120 seconds 

TIME STEP : 30 seconds 


ELAPSED TIME 00:30 
TWO VIEW APPROACH CENTER LINE APPROACH 


MAX. HEIGHT 3; 70 46 METERS 

MAX. DISTANCE 3; 48 49 METERS 

VOLUME OF SMOKE ; 41338 73638 CU. METERS 

PRODUCTION RATE : 82676 147277 CU. METERS PER MINUTE 


PLUME TRAJECTORY 


TEST NUMBER : 5 
ELAPSED TIME 00:30 
COMPASS HEADING 0 degrees 


DISTANCE (M) HEIGHT (M) 
0 0 
9.936396 8.889222 
19.87281 17.9468 
29.8092 27.17275 
39.74561 36.56705 
49.68201 46.12972 
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ELAPSED TIME 01:00 


TWO VIEW APPROACH CENTER LINE APPROACH 
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MAX. HEIGHT : 130 114 ME 
MAX. DISTANCE : 122 191 eens 
VOLUME OF SMOKE : 475385 541395 CU. METERS 
PRODUCTION RATE : 868095 935512 CU. METERS PER MINUTE 
PLUME TRAJECTORY 
TEST NUMBER : 5 
ELAPSED TIME 01:00 
COMPASS HEADING 0 degrees 
DISTANCE (M) HEIGHT (M) 
0 0 
11.17844 8.899744 
22.35688 17.98889 
33.53535 27.26745 
44.71379 36.73541 
55.89225 46.39278 
67.07063 56.23949 
91.91166 67.0859 
116.7527 78.35321 
141.5937 90.04142 
166.4347 102.1505 
191.2757 114.6806 
ELASPED TIME 01:30 
TWO VIEW APPROACH CENTER LINE APPROACH 
MAX. HEIGHT : 200 172 METERS 
MAX. DISTANCE : 198 308 METERS 
VOLUME OF SMOKE : 1974481 2470698 CU. METERS 
PRODUCTION RATE : 2998192 3858606 CU. METERS PER MINUTE 
PLUME TRAJECTORY 
TEST NUMBER : 5 
ELAPSED TIME 01:30 
COMPASS HEADING 0 degrees 
DISTANCE (M) HEIGHT (M) 
1) i) 
9.93641 7.619325 
19.87276 15.38296 
29.80917 23.2909 
39.74558 31.34315 
49.68193 39.53971 
59.61834 47.88058 
59.6184 47.88058 
109.3004 69.50726 
158.9824 92.81756 
208.6644 117.8115 
258.3464 144.489 
308.0285 172.8501 


ELAPSED TIME 02:00 
TWO VIEW APPROACH CENTER LINE APPROACH 


MAX. HEIGHT 3; 310 302 METERS 

MAX. DISTANCE ; 350 541 METERS 

VOLUME OF SMOKE ; 6536319 7025034 CU. METERS 

PRODUCTION RATE 3; 9123675 9108671 CU. METERS PER MINUTE 


e 


PLUME TRAJECTORY 


TEST NUMBER : 5 
ELAPSED TIME 02:00 
COMPASS HEADING O degrees 


DISTANCE (M) HEIGHT (M) 
0 0 
9.93641 10.15912 
19.87282 20.51065 
29.80923 31.0546 
39.74564 41.79095 
49.68205 52.71972 
59.61834 63.84078 
119.2367 81.3259 
178.8551 100.2541 
238.4735 120.6254 
298.0919 142.4398 
357.7103 165.6973 
417.3284 190.3978 
442.1694 211.5748 
467.0104 233.3532 
491.8514 255.7328 
516.6924 278.7136 
541.5333 302.2958 
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agreement was obtained between corresponding 35 mm photographs and digitized images used to 
estimate plume volume and trajectory. In Test 7, at 600 s after ignition, the leading edge 
of the smoke plume had risen to 750 m and traveled 1800 m downwind from the pag. .The 
estimated total volume of the smoke plume at 600 s after ignition was 3.8 x 10¢ m,. At that 
time, the total volume of the smoke plume was increasing at a rate of 3.2 x 10 m/s. Limited 
evaluation indicates that the MS-DOS based method provides reasonably accurate estimates of 
visible smoke plume geometry in the near-field. Field accuracy depends on plume size, wind 
speed and direction, and the resolution of the equipment. 
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